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FOREWORD 


Dipole  Pride  26  test  program  participants  included  support  elements  from  the 
Department  of  Energy  (DOE),  the  National  Oceanographic  and  Atmospheric  Administration 
Air  Resources  Laboratory  (NOAA  ARL),  the  Department  of  the  Army  (DA),  and  several 
contractors.  Test  program  oversight  and  logistical  support  were  provided  by  the  DOE 
Nevada  Operations  Office.  The  NOAA  ARL  Special  Operations  and  Research  Division 
(ARL/SORD)  provided  meteorological  support,  while  the  NOAA  ARL  Field  Research  Division 
(ARL/FRD)  performed  the  tracer  gas  measurements  using  sulfur  hexafluoride  detectors. 

The  U.S.  Army  Dugway  Proving  Ground  West  Desert  Test  Center  (DPGA/VDTC)  designed 
and  managed  the  test  program,  operated  the  puff  dissemination  system,  and  provided 
micrometeorological  measurements.  Test  program  contractors  included  Logicon  RDA,  The 
Aerospace  Corporation,  and  Bechtel  Nevada.  Logicon  contributed  technical  expertise  to 
various  aspects  of  test  design.  The  Aerospace  Corporation  operated  remote  sensing 
equipment  to  track  the  puffs,  and  Bechtel  Nevada  provided  onsite  logistical  support. 


Blank  Page 


ii 


ACKNOWLEDGMENTS 

The  Phase  II  Dipole  Pride  26  (DP26)  Defense  Special  Weapons  Agency  (DSWA)  model 
validation  program  included  support  elements  from  the  Department  of  Energy  (DOE)  Nevada 
Test  Site  staff  managed  by  MAJ  Paul  Loomis.  Test  program  participants  included  the 
National  Oceanographic  and  Atmospheric  Administra-tion  Air  Resources  Laboratory  (NOAA 
ARL),  the  Department  of  the  Army  (DA),  DSWA,  and  several  contractors.  Dr.  Darryl 
Randerson,  Director  of  the  NOAA  ARL  Special  Operations  and  Research  Division  (SORD) 
provided  operating  space  at  the  Yucca  Flat  Weather  Station  and  laboratory  space  for  the 
gas  laboratory  in  Mercury,  Nevada.  Mr.  Ray  Dennis,  Senior  SORD  Meteorological  techni¬ 
cian  and  his  crew  provided  invaluable  day-to-day  meteorological  support.  Drs.  David 
George  and  Thomas  Watson  of  the  NOAA  ARL  Field  Research  Division  (FRD)  managed  the 
field  crews  and  gas  laboratory  personnel  that  collected  the  sulfur  hexafluoride  sampler 
data.  LTC  A.J.  Kuehn  represented  DSWA  at  the  test  site,  with  Drs.  Gary  Ganong  and 
William  Espander  from  Logicon  RDA  providing  valuable  technical  assistance.  Puff  imagery 
and  FTIR  spectronomy  were  provided  by  The  Aerospace  Corporation  under  the  direction  of 
Dr.  Kenneth  Herr.  Charles  Birdsong,  Christopher  Woldruff,  and  David  Petrie  of  the  DPG 
West  Desert  Test  Center  (WDTC)  Test  Operations  Division  operated  the  dissemination 
system,  provided  puff  source  dimensions,  and  collected  micrometeorological  data.  Drs. 
Tom  Watson,  Steve  Hanna,  and  Darryl  Randerson  provided  external  review  of  this  report. 
Mrs.  Susan  Gross  of  the  WDTC  Meteorology  &  Obscurants  Division  provided  word 
processing  support. 


iii 


Blank  Page 


EXECUTIVE  SUMMARY 


Phase  II  of  the  Defense  Special  Weapons  Agency  (DSWA)  Transport  and  Dispersion 
Model  Validation  Program  (Dipole  Pride  26)  was  conducted  at  Yucca  Flat  on  the  Nevada 
Test  Site  in  November  1996.  The  Phase  II  test  objective  was  to  acquire  a  data  base  for  the 
validation  of  integrated  mesoscale  wind  field  and  dispersion  models,  in  particular  the 
Hazard  Prediction  and  Assessment  Capability  (HPAC)  model  suite.  This  objective  was 
achieved  by  releasing  tracer  gas  (sulfur  hexafluoride)  puffs,  with  downwind  tracer  sampling 
at  distances  ranging  to  20  km  along  with  extensive  measurements  to  document  meteoro¬ 
logical  conditions.  DSWA  sponsored  the  Phase  II  niodel  validation  test  series  as  part  of  its 
counterproliferation  model  development  and  validation  effort. 

A  total  of  23  puff  releases  were  completed  during  the  trials  series,  including  5  into 
nocturnal  drainage  flows.  Significant  terrain  effects  observed  during  puff  travel  through 
Yucca  Flat  included  pooling  at  Yucca  Lake  during  drainage  flows,  and  upslope  (katabatic) 
flows  as  the  mountains  on  either  side  of  Yucca  Flat  received  solar  heating.  Crosswind  puff 
dispersion  was  sampled  at  three  sampling  lines  located  along  roads  2  to  20  km  downrange 
of  the  source.  Alongwind  puff  growth  was  also  measured  along  the  middle  sampling  line, 
but  interpretation  of  these  measurements  is  complicated  by  disseminator  system  leaks. 

This  report  presents  lateral  (Oy)  and  alongwind  (a,)  puff  dispersion  summaries  ob¬ 
tained  from  tracer  concentration  measurements,  plus  detailed  explanations  of  puff 
dissemination  and  sigma  calculation  procedures.  Remote  imagery  was  used  to  obtain 
information  on  vertical  dispersion  and  the  puff  centroid  position.  Results  from  the  fixed 
sampling  lines,  which  include  whole  air  sampler  data  and  continuous  analyzer  data,  are 
available  on  a  compact  disk  along  with  the  supporting  meteorological  and  micrometeor- 
ological  data.  Report  documentation  includes  the  puff  dissemination  conditions,  downwind 
dispersion,  atmospheric  stability,  site  roughness,  and  micrometeorological  summaries 
needed  for  dispersion  model  validation. 
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SECTION  1 ,  INTRODUCTION 


Responding  to  a  shortfall  in  downwind  hazard  effects  modeling  identified  during 
Operation  Desert  Storm,  the  Defense  Special  Weapons  Agency  (DSWA)  has  developed  and 
is  evaluating  the  performance  of  the  Hazard  Prediction  and  Assessment  Capability  (HPAC) 
model  suite,  HPAC  includes  a  diagnostic  wind  field  model  coupled  with  an  atmospheric 
transport  and  dispersion  model  to  provide  decision  makers  with  a  tool  for  predicting 
windborne  hazards  that  arise  as  a  consequence  of  toxic  materials  released  into  the 
atmosphere.  The  HPAC  development  effort  began  with  the  recognition  that  the  ensemble 
mean  predictions  available  from  the  current  generation  of  atmospheric  dispersion  models 
used  for  chemical  and  biological  (CB)  hazard  assessments  are  of  limited  operational  utility. 
Ensemble  mean  predictions  provide  an  estimate  of  the  concentrations  or  dosages  that 
would  be  obtained  by  averaging  the  results  for  a  large  number  of  replications  of  the  same 
release,  but  include  no  information  about  the  distribution  of  single  event  concentrations  or 
dosages.  Current  hazard  assessment  models  also  cannot  estimate  the  probability  of 
exceeding  critical  hazard  thresholds.  These  deficiencies  complicate  the  interpretation  of  CB 
hazard  model  predictions  and  limit  their  usefulness  in  operational  military  situations. 

Given  the  limited  operational  usefulness  of  ensemble  mean  predictions,  DSWA 
identified  the  need  for  a  probabilistic  atmospheric  dispersion  model.  In  addition  to  providing 
a  prediction  of  ensemble  mean  concentrations  or  dosages,  a  probabilistic  dispersion  model 
uses  higher  order  statistical  terms  to  define  the  probability  that  dosages  or  peak  concentra¬ 
tions  at  points  of  interest  exceed  some  critical  value  specified  by  the  user.  This  probabilis¬ 
tic  output  requires  that  the  model  be  able  to  predict  the  concentration  or  dosage  cumulative 
distribution  function  (CDF),  which  in  turn  requires  that  the  model  be  able  to  predict  both 
the  concentration  (or  dosage)  means  and  variances.  The  Second  Order  Closure  Integrated 
Puff  (SCIPUFF®)  model  (Titan  Corporation,  1996),  which  predicts  the  required  means  and 
variances  through  a  second-order  closure  solution  of  the  advection-diffusion  equation,  is 
the  HPAC  dispersion  model  component. 

SCIPUFF  uses  a  generalized  Gaussian  tensor  to  describe  puff  concentrations.  The 
model's  derivation  begins  by  integrating  the  conservation  equations  to  obtain  differential 
equations  for  the  concentration  moments  that  explicitly  include  wind-shear  effects. 
Second-order  closure  is  used  to  relate  these  higher-order  terms  to  turbulence  parameters 
such  as  velocity  cross-correlations  and  turbulence  length  scales.  In  contrast  to  current 
generation  Gaussian  puff  models,  SCIPUFF’s  second-order  closure  methodology  yields 
theoretically  consistent  predictions  of  dosage  and  concentration  variances  in  addition  to 
their  ensemble  means.  These  means  and  variances  are  applied  in  a  clipped-normal 
probability  distribution  to  predict  concentration  and  dosage  CDFs. 

In  the  fall  of  1 995,  DSWA  contacted  the  Dugway  Proving  Ground  (DPG)  West 
Desert  Test  Center  (WDTC)  Meteorology  &  Obscurants  Division  (WD-M)  for  assistance  in 
the  development  and  implementation  of  a  transport  and  dispersion  model  validation 
program.  WD-M  convened  a  meeting  of  atmospheric  transport  and  dispersion  modeling 
experts  in  February  1 996  to  recommend  a  model  validation  program.  Attendees  included 
representatives  from  the  Joint  Services,  the  National  Oceanic  and  Atmospheric  Administra¬ 
tion  (NOAA),  and  the  trinational  (U.S.,  UK,  and  Canada)  Technical  Panel  9  of  The  Technical 
Cooperation  Program  (TTCP)  Subgroup  E  on  CB  Defense.  Recognizing  the  limitations  of  the 
data  from  CB  weapons  tests  and  previous  field  dispersion  experiments,  the  attendees 
agreed  that  new  data  sets  and  model  validation  procedures  were  needed  to  validate 
SCIPUFF's  probabilistic  output.  Consequently,  the  attendees  suggested  conceptual  designs 
for  a  series  of  field  tests  to  create  high  resolution  puff  data  sets  accompanied  by  detailed 
meteorological  documentation. 
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The  DSWA  Model  Validation  Program  has  had  two  phases  to  date.  Phase  I, 
conducted  9-26  September  1 996  at  DPG,  consisted  of  observations  of  short-range  (200  to 
1200  m)  puff  releases  for  CDF  evaluations.  Phase  II  considered  transport  and  diffusion  to 
mesoscale  distances  (10  to  20  km).  The  Phase  II  subtest,  known  as  Dipole  Pride  26 
(DP26),  consisted  of  a  series  of  puff  releases  conducted  at  the  Nevada  Test  Site  on  4-21 
November  1996.  This  report  documents  the  Phase  II  (DP26)  test  program  and  its  results. 
Detailed  sampler  data  and  supporting  meteorological  and  micrometeorological  measure¬ 
ments  are  also  available  on  magnetic  and/or  optical  media. 
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SECTION  2.  TEST  DESCRIPTION 


2.1  TEST  SITE 

Yucca  Flat  (37°  N,  116°  W)  is  a  north-south  oriented  basin  30  km  in  length  and  1 2 
km  in  width  surrounded  on  all  sides  by  higher  terrain.  Yucca  Lake,  at  a  height  of  1 195  m 
above  mean  sea  level  (MSL),  is  a  seasonally  dry  lake  bed  at  the  south  end  of  the  basin  that 
forms  the  lowest  part  of  the  Yucca  Flat.  The  basin  is  surrounded  by  mountains  extending 
to  1 800  m  MSL  or  higher  on  the  west  through  north  and  to  1 500  m  MSL  or  higher  on  the 
east.  The  basin  slopes  upward  toward  the  north  at  an  average  angle  of  0.3°.  Wind  flow 
through  Yucca  Flat  is  strongly  influenced  by  the  surrounding  terrain.  Passes  to  the 
northeast  and  south  form  the  main  ventilation  channels  through  Yucca  Flat.  Early  morning 
drainage  flows  from  the  north  cause  pools  of  cold  air  to  accumulate  over  Yucca  Lake. 

Solar  heating  of  higher  terrain  at  the  north  end  of  Yucca  Flat  generally  draws  a  southerly 
flow  through  the  basin  during  daylight  hours. 

The  principal  test-related  facilities  within  Yucca  Flat  included  the  Yucca  Flat 
Weather  Station  (designated  UCC  for  meteorological  reports)  on  the  western  edge  of  Yucca 
Lake;  the  Buster-Jangle  Yankee  intersection  (designated  BJY),  a  position  near  the  intersec¬ 
tion  of  Rainer  Mesa  Road  and  Mercury  Highway;  and  a  network  of  meteorological  data 
(MEDA)  stations  within  and  around  Yucca  Flat.  These  locations  are  shown  on  Figure  1 , 
and  their  latitude  and  longitude  positions  are  given  in  Table  1.  The  ARL/SORD  meteoro¬ 
logical  building  at  UCC  served  as  both  the  test  program  command  post  and  a  release  point 
for  radiosonde  and  pilot  balloon  (pibal)  flights.  A  pibal  station  and  wind  profiling  radar  were 
also  stationed  at  BJY.  Tower-mounted  sonic  anemometer/thermometers  (sonics)  were 
operated  near  BJY  and  at  a  site,  designated  YFW  in  this  report,  that  is  1 00  m  east  of  UCC. 


Table  1.  Dipole  Pride  26  Test  Site  and  MEDA  Station  Locations. 


Position 

Latitude  (°N) 

Longitude  (®W) 

Area  1  (MEDA  1 ,  elevation  1 265  m) 

37.0275 

116.0917 

Area  2  (MEDA  2,  elevation  1341  m) 

37.1392 

116.1058 

Area  3  (MEDA  3,  elevation  1 207  m) 

37.0042 

116.0317 

Area  9  (MEDA  9,  elevation  1 290  m) 

37.1358 

116.0400 

BJY  (MEDA  1 7,  elevation  1 244  m) 

37.0625 

116.0525 

DAF  (MEDA  28,  elevation  1 107  m) 

36.8925 

116.0375 

YFW  (MEDA  6,  elevation  1 1 95  m) 

36.9583 

116.0467 

MON  (MEDA  10,  elevation  1570  m) 

36.9400 

116.0792 

CSE  (Cane  Springs  Road  East) 

36.8512 

115.9985 

N2  (North  Dissemination  Site  2) 

37.1586 

1 16.0967 

N3  (North  Dissemination  Site  3) 

37.1500 

116.0625 

S2  (South  Dissemination  Site  2) 

36.9570 

116.0498 

S3  (South  Dissemination  Site  3) 

36.9512 

116.0100 

Sampler  Site  101  (West  End  North  Line) 

37.1383 

116.1230 

Sampler  Site  130  (East  End  North  Line) 

37.1227 

1 16.0410 

Sampler  Site  201  (West  End  Middle  Line) 

37.0557 

116.0920 

Sampler  Site  230  (East  End  Middle  Line) 

37.0493 

1 16.0090 

Sampler  Site  301  (West  End  South  Line) 

36.9906 

116.0930 

Sampler  site  330  (East  End  South  Line) 

36.9954 

116.0160 
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2.2  FIELD  TEST  DESIGN 


The  Dipole  Pride  26  (DP26)  field  test  was  designed  to  use  fixed  sampling  lines  and  a 
mobile  dissemination  system  to  document  puff  dispersion  utilizing  the  predominant  north- 
south  flow  through  Yucca  Flat.  Six  dissemination  positions,  three  on  the  north  side  and 
three  on  the  south  side  of  Yucca  Flat,  were  selected  for  the  truck-mounted  disseminator  to 
accommodate  wind  direction  variations.  Only  sites  S2,  S3,  N2,  and  N3  were  actually  used 
for  dissemination. 

The  dissemination  system  was  built  by  Consumers  Pipe  in  Las  Vegas,  NV.  It 
consisted  of  gas  cylinders,  the  tracer  gas  sulfur  hexafluoride  .(SFg),  a  compressor,  actua¬ 
tors,  fill  hoses,  and  plumbing  mounted  on  the  bed  of  a  government-furnished  5-ton  truck. 
Gas  was  released  through  two  cylinders  with  fast-acting  solenoid-actuated  butterfly  valves. 
Each  release  was  completed  within  2  seconds  of  valve  opening.  The  cylinders  were 
designed  to  be  operated  either  singly  or  together.  Each  cylinder  had  a  volume  of  0.1 5  m^ 
and  held  approximately  10  kg  of  SFg  when  filled  to  the  design  pressure  of  150  lbs  per 
square  inch  (105,465  kg/m^).  Table  2  gives  the  calculated  quantities  of  SFg  released  for 
each  trial,  and  Appendix  B  describes  the  methodology  used  to  calculate  release  quantities. 

The  released  SFe  puffs  were  sampled  using  six  TGA-4000  Tracer  Gas  Analyzers  and 
ninety  whole  air  samplers  dispersed  along  three  sampling  lines.  Technicians  operated  the 
TGA-4000S,  which  were  installed  within  half-ton  cargo  vans.  These  technicians  also 
serviced  whole  air  samplers  mounted  on  posts  1.5  m  above  ground  level  (AGL)  along 
sampling  line  roads.  The  fast-response  TGA-4000s  provided  real  time  ihformation  on  the 
arrival  and  departure  of  each  puff  and  measured  the  high-frequency  (4-Hz)  variations  of  the 
gas  concentration  field  within  the  passing  puffs.  The  whole  air  samplers  provided  only 
coarse  (15-min)  time  resolution,  but  the  time-averaged  gas  concentration  measurements 
from  the  30  whole  air  samplers  stationed  along  east-west  roads  crossing  Yucca  Flat 
provided  the  spatial  resolution  required  for  puff  lateral  dimension  calculations.  Additional 
whole  air  samplers  were  positioned  away  from  the  planned  puff  trajectory  as  background 
SFg  references,  and  dual  samplers  were  mounted  at  some  sampling  stations  for  quality 
control. 

DP26  sampling  lines  were  established  along  the  following  roads:  Road  2-04 
between  Orange  Road  and  Mercury  Highway  (Line  1,  sampler  positions  101  through  130); 
Pahute  Mesa  Road  between  Orange  Road  and  Orange  Blossom  Road  (Line  2,  sampler 
positions  201  through  230);  Mine  Mountain  Road  between  Tippipah  Highway  and  Orange 
Blossom  Road  (Line  3,  sampler  positions  301  through  330).  The  TGA-4000s  were 
stationed  at  1 500-m  intervals  along  Sampling  Line  2  (Pahute  Mesa  Road)  to  maximize  the 
possibility  of  intercepting  portions  of  the  cloud  crossing  this  line.  Mean  whole  air  sampler 
spacing  along  Lines  1 ,  2,  and  3  were  259,  256,  and  237  m,  respectively.  Sampler 
spacing,  line  length,  and  number  was  chosen  as  a  compromise  between  four  constraints: 

(1)  containing  the  entire  puff  width  within  a  sampling  line;  (2)  intercepting  the  puff  with  at 
least  six  samplers;  (3)  remaining  within  the  confines  of  existing  roads;  and  (4)  not 
exceeding  technician  servicing  capabilities.  This  compromise  optimized  sampling  along  Line 
2,  which  was  roughly  10  km  from  the  north  and  south  dissemination  positions.  Pre-trial 
modeling  suggested  a  sampler  spacing  of  250-300  m  as  the  best  compromise  between  the 
first  and  second  constraints.  The  third  and  fourth  constraints  limited  sampler  line  length 
and  the  number  of  positions  within  a  sampling  line. 

Puff  infrared  (IR)  imagery  was  obtained  by  WDTC  personnel  at  the  time  of  release, 
and  the  dispersing  puff  was  tracked  downrange  using  an  array  of  mobile  imagers  operated 
by  The  Aerospace  Corporation.  The  principal  purposes  of  this  imagery  were  to  document 
initial  puff  dimensions  and  to  resolve  in  time  and  three-dimensional  space  the  puff  centroid 
positions  as  these  centroids  crossed  the  sampling  lines.  Operators  of  the  Aerospace 
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imagers  selected  their  positions  prior  to  puff  release  based  on  release  location  and  wind 
conditions. 

The  Dipole  Pride  26  test  series  was  designed  to  complement  puff  dispersion 
measurements  with  extensive  meteorological  documentation.  This  documentation  included 
surface-based  and  upper  air  measurements  and  weather  analysis  and  forecasting  support 
from  the  ARL/SORD  facility  in  Las  Vegas,  NV.  Surface  measurements  included  wind, 
temperature,  humidity,  and  pressure  data  obtained  over  1 5-min  averaging  periods  from  an 
array  of  MEDA  stations  located  in  and  around  Yucca  Flat.  These  measurements  were 
supplemented  with  temperature  and  momentum  flux  and  turbulence  information  obtained 
using  sonic  anemometer/thermometers  mounted  on  towers  at  BJY  and  YFW.  Pibal  and 
radiosonde  flights  provided  wind  and  thermodynamic  profiles  through  the  lower  tro¬ 
posphere.  ARL/SORD  reactivated  its  radiosonde  site  UUC  at  the  Yucca  Flat  Weather 
Station  during  DP26  to  launch  radiosonde  flights  every  3  hours  while  trials  were  in  process. 
Pibals  were  also  launched  every  hour  from  BJY  and  UUC  (except  during  radiosonde 
flights),  and  occasionally  from  Cane  Springs  Road.  A  924-MHz  wind  profiling  radar 
positioned  near  BJY  provided  intermittent  wind  profiles  during  testing.  MEDA  station, 
pibal,  radiosonde,  and  profiler  data  are  available  on  floppy  disks  or  CD  ROM.  The  sonic 
anemometer/thermometer  measurements  require  extensive  processing  to  produce  useful 
statistical  results.  These  statistical  results  are  presented  in  this  report,  while  the  basic  data 
are  available  from  DPG  by  special  request. 

2.3  TEST  INSTRUMENTATION 

2.3.1  TGA-4000  Continuous  SFc  Analyzer 

The  TGA-4000  is  designed  to  provide  real-time  fast  response  (better  than  1  Hz)  SFg 
concentration  measurements.  It  consists  of  a  catalytic  reactor,  a  dryer,  and  an  electron 
capture  detector  (ECD)  to  measure  halogenated  compounds  such  as  SFg.  Because  the 
detector  responds  to  oxygen,  the  TGA-4000  uses  the  catalytic  reactor  to  convert  molecular 
oxygen  to  water  and  a  dryer  subsequently  removes  the  water  from  the  airstream.  Any  SFe 
above  a  picogram  threshold  remaining  in  the  airstream  after  these  operations  produces  a 
voltage  proportional  to  its  concentration  as  the  airstream  passes  through  the  ECD.  The 
TGA-4000,  its  associated  plumbing,  and  a  data  acquisition  and  display  system  fit  into  a  ’72- 
ton  cargo  van,  which  can  be  driven  into  sampling  position.  The  van  also  serves  as  a  power 
source  for  the  TGA-4000.  Scientech,  Inc.  of  Pullman,  WA,  who  built  the  TGA-4000, 
report  a  bench-top  noise  level  of  5  part  per  trillion  by  volume  (pptv)  and  a  response  time  of 
0.86  s  (Benner  and  Lamb,  1985).  A  4-Hz  data  acquisition  rate  was  used  with  the  TGA- 
4000  instruments  during  DP26.  Further  details  on  the  TGA-4000  can  be  found  in  Watson 
et  al.  (1998)  or  Bowers  et  al.  (1994). 

2.3.2  Whole  Air  Samplers 

Whole  air  samplers,  also  known  as  sequential  bag  samplers  (Bowers  et  al,  1994), 
consist  of  a  cartridge  containing  a  programmable  microprocessor  and  twelve  1-1  Tedlar™ 
bags,  with  each  bag  connected  to  a  small  air  pump.  A  "D"  cell  battery  provides  power  for 
the  whole  air  sampler  unit.  The  cartridge  is  mounted  within  a  waxed  cardboard  box  that 
fits  onto  a  mounting  bracket.  During  DP26  each  sampler  was  mounted  on  a  post  at  a 
height  of  1.5  m  above  ground  level  (AGL).  The  sampler  can  be  programmed  to  sequentially 
fill  each  of  the  bags  over  time  periods  ranging  from  10  min  to  several  hours,  thereby 
providing  a  concentration  measurement  integrated  over  the  selected  sampling  period.  A 
1 5-min  sampling  period  was  used  during  DP26  as  the  best  compromise  between  the  desire 
for  maximum  time  resolution  and  the  need  to  sample  during  the  3-hour  duration  of  a  typical 
test.  The  bag  1  fill  sequence  was  programmed  to  begin  at  puff  dissemination  time  on  the 
two  sampling  lines  closest  to  the  disseminator,  and  was  delayed  for  30  min  on  the  furthest 
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sampling  line.  Once  the  fill  sequence  on  bag  1 2  was  completed,  the  cartridges  containing 
sealed  bags  were  dismounted  and  taken  to  a  field  gas  laboratory  set  up  in  Mercury,  NV  for 
separation  by  species.  Chemical  species  separation  was  performed  using  a  gas  chroma¬ 
tograph  and  an  electron  capture  detector.  Further  details  on  whole  air  sampler  operation 
during  DP26  are  given  by  Watson  et  al  (1998). 

2.3.3  Sonic  Anemometer/Thermometers 

A  sonic  anemometer/thermometer  (sonic)  consists  of  a  transducer  array  containing 
paired  sets  of  ultrasonic  transmitters  and  receivers,  a  system  clock,  and  circuitry  designed 
to  measure  intervals  of  time  from  the  transmission  to  the  reception  of  sound  pulses 
traveling  between  transducer  pairs.  The  sonics  used  during  Dipole  Pride  26  were  the 
Applied  Technologies,  Inc.  3-axis  ATI  Model  RSWS-201/3A,  which  provides  three  dimen¬ 
sional  (u,  V,  and  w)  wind  components  and  speed  of  sound.  These  sonics  were  rnounted  at 
10  and  14  m  AGL  on  towers  near  the  Buster-Jangle  Yankee  (BJY)  intersection  and  at  a 
tower  designated  YFW  located  100  m  east  of  the  Yucca  Flat  Weather  Station  (UCC). 

Sonic  locations  were  confined  to  these  sites  by  the  need  for  a  suitable  mounting  tower  and 
1 1 0  VAC  power.  The  3-axis  sonics  resolve  all  components  of  the  three-dimensional  wind 
vector.  The  speed  of  sound  measurements  were  converted  to  sonic  temperature,  which  is 
essentially  equivalent  to  the  virtual  temperature.  With  a  data  rate  of  1 0  Hz  and  an  acoustic 
pathlength  between  transducers  of  1 5  cm,  the  sonics  provide  sufficient  temporal  and 
Spatial  resolution  to  measure  mean  wind  plus  the  fluctuating  components  needed  to  define 
turbulence  intensities  and  the  fluxes  of  heat  and  momentum.  Procedures  described 
American  Society  for  Testing  of  Materials  (ASTM)  standard  practice  for  obtaining  wind  and 
temperature  measurements  from  sonics  (ASTM  1997A)  were  used  to  obtain  wind  compo¬ 
nent  and  speed  of  sound  measurements  to  within  ±3  cm/s. 


2.3.4  MEDA  Stations 

ARL/SORD  operates  a  continuousiy  monitoring  network  of  remote  meteorological 
stations  (MEDA  stations)  across  the  Nevada  Test  Site.  Each  station  in  the  MEDA  network 
reports  a  measurement  of  temperature,  pressure,  humidity,  wind  speed,  and  wind  direction 
every  1 5  min.  There  is  a  time  tag  at  the  end  of  each  MEDA  data  block.  For  example,  dap 
collected  during  the  period  0230  -  0245  Pacific  Standard  Time  (PST)  are  reported  at  0245 
PST  With  the  exception  of  wind  sensors  mounted  at  10  m,  MEDA  station  instruments  are 
mounted  at  2  m  AGL.  MEDA  station  data  are  transmitted  via  radio  to  a  centralized  data 
collection  point.  The  reported  temperature,  pressure,  and  humidity  are  the  last  availaWe 
readings  during  each  1 5-min  reporting  interval.  Wind  speed  and  direction  are  obtained  over 
a  5-min  averaging  period  immediately  prior  to  the  reporting  time,  with  maximum  and 
minimum  speeds  taken  from  1-s  readings  during  the  entire  1 5-min  period.  The  MEDA 
station  data  displays  in  the  Yucca  Flat  Weather  Station  served  as  the  primary  source  of 
wind  information  for  DP26  test  conduct. 


2.3.5  Radiosonde 

A  radiosonde  system  consists  of  a  balloon-borne  instrument  package  that  rises 
through  the  atmosphere,  providing  profiles  of  wind,  temperature,  humidity,  and  height  at 
10-s  intervals  throughout  the  flight.  Data  packets  from  this  instrument  are  transmitpd  via 
radio  link  to  a  base  station  where  they  are  logged  and  subjected  to  quality  control.  “Ch 
profile  is  reduced  to  measurements  at  standard  pressure  levels  (the  mandatory  levds)  and 
significant  inflection  points  (significant  levels).  The  radiosonde  system  used  at  N '  S  |s  the 
automatic  radio-theodolite  (ART),  which  tracks  the  balloons  with  a  ground-based  radio 
theodolite.  The  ART  tracks  the  ascending  balloon,  providing  elevation  and  azimuth  angle 
readings  which,  after  time-synchronization  with  pressure  readings,  are  converted  to  profiles 
of  wind  speed  and  direction.  The  ART  system  used  during  DP26  was  located  at  the  Yucca 
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Flat  Weather  Station  (UCC).  It  provided  wind  and  thermodynamic  profiles  every  3  hours 
during  the  test  program.  The  ART  measurement  uncertainties  are  ±0.6  °C  for  temperature, 
±10%  for  relative  humidity,  and  ±1  m/s  for  winds  (RCC-MG,  1992). 

2.3.6  Pibal 


The  pilot  balloon  (pibal)  is  an  optically-tracked  free  balloon  used  to  obtain  profiles  of 
wind  speed  and  direction.  Hundred  gram  (100-g)  pibals  provided  boundary  layer  wind 
profiles  during  this  test  program.  When  filled  to  its  design  lift  weight,  a  100-g  pibal  has  an 
ascent  rate  that  is  large  in  comparison  with  typical  atmospheric  vertical  motions.  Standard 
tables  are  used  to  relate  a  pibal’s  flight  time  to  its  height  AGL.  Optical  tracking  with  a 
theodolite  provides  azimuth  and  elevation  readings  taken  at  30-s  intervals.  These  readings, 
combined  with  tabulated  height  versus  time  data,  provide  sufficient  information  to  calculate 
layer-averaged  wind  speeds  and  directions.  Pibal  wind  profiles  are  typically  accurate  to 
within  ±2  m/s  (RCC-MG,  1992).  The  Digital  Pibal  (DIGIPI)  Systems  used  at  NTS  feature 
shaft  encoders  to  digitally  record  theodolite  angles  every  30  s.  The  angular  data  are  stored 
in  a  microcomputer  linked  to  a  communications  device  that  transmits  these  data  to  the 
central  computer  during  a  polling  sequence.  DIGIPI  units  were  stationed  at  BJY  and  UCC, 
and  occasionally  at  Cane  Springs  Road  (CSE)  in  Frenchman  Flat,  providing  hourly  wind 
profiles  during  the  test  program. 

2.3.7  Infrared  Imaging  Radiometers 

Infrared  imaging  radiometers  (imagers)  are  passive  optical  devices  sensitive  to  IR 
energy  in  the  8-  to  1 2-mm  portion  of  the  IR  spectrum.  Because  SFg  has  a  distinct 
absorption  band  in  this  portion  of  the  spectrum,  the  passage  of  a  SFg  puff  across  the 
imager's  field  of  view  registers  as  a  temperature  change  when  compared  with  the  back¬ 
ground  image.  IR  imagers  were  used  to  determine  initial  puff  dimensions  (source  size)  and 
to  track  the  puffs  as  they  traversed  the  sampling  lines.  Source  size  characterization  was 
done  using  an  Inframetrics  Model  600L  Imaging  Radiometer,  which  has  a  typical  thermal 
sensitivity  of  0.05  “C,  a  scan  rate  of  50  Hz,  and  a  7-bit  (128  levels)  image  resolution. 

Three  Sterling-cooled  Agema  Thermovision  900  digital  IR  imagers,  each  equipped  with  a 
narrowband  filter  centered  around  the  main  SFg  absorption  band,  were  used  to  monitor  SFg 
puff  travel  across  the  sampling  lines.  These  imagers  are  characterized  by  a  1 5-Hz  scan 
rate  (non-interlaced)  and  a  thermal  sensitivity  of  0.08  X.  The  spectral  resolution  and  field 
of  view  varied  with  the  filters  and  lenses  used.  These  images  captured  information  on  puff 
concentration  column  density  and  centroid  height. 

2.3.8  Fourier  Transform  Infrared  Spectrometer 

The  Fourier  transform  infrared  spectrometer  (FTIR)  is  a  remote  imaging  device 
designed  to  measure  the  scene  spectrum  within  the  8-  to  1 2-mm  band.  Because  all 
materials  that  absorb  IR  energy  have  unique  absorption  band  signatures,  an  FTIR  is  able  to 
resolve  the  contents  of  a  dispersing  cloud  if  a  sufficient  signature  is  present.  These 
instruments  are  most  useful  when  the  background  within  their  fields  of  view  is  constant 
and  the  target  tracer  cloud  is  easily  distinguishable  from  the  background  thermal  signature. 
Consequently,  the  FTIRs  were  oriented  toward  the  (relatively  cold)  clear  sky  to  detect  the 
arrival  and  passage  across  the  field  of  view  of  the  (relatively  warm)  SFg  puffs.  Sulfur 
hexafluoride  is  an  ideal  tracer  because  it  exhibits  a  sharp  peak  at  an  inverse  wavelength  of 
950  cm'\  When  processed  through  a  "special  ratio"  algorithm  (Polak  et  al.,  1995) 
designed  to  remove  background  and  noise,  the  SFg  transmission  or  absorption  peak 
emerges  and  can  be  used  for  quantitative  estimation  of  puff  column  density. 
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The  FTIR  supporting  DP26  was  an  Intillitec  M21  chemical  agent  detector.  This 
instrument  features  a  spectral  resolution  of  1.5  cm  \  a  5.25-Hz  scan  rate,  a  25-milliradian 
field  of  view,  and  a  sensitivity  of  1 .5  x  10  ®  W  cm  ^  sr  Vcm  ’.  Spectrometers  were 
mounted  in  the  Aerospace  Corporation  Ram  Van  which  was  positioned  to  intercept  the 
passing  SFe  puff  centroid,  and  in  the  Aerospace  Tonka  Van,  which  followed  the  puff  and 
made  traverses  through  it  as  it  dispersed  downrange.  Reports  from  these  two  vans 
provided  real-time  information  on  puff  location  and  dispersion,  and  follow-on  analysis  of 
these  data  could  provide  valuable  puff  position  information. 

2.3.9  Dissemination  System 

The  SFg  gas  release  system  consisted  of  two  vertically-mounted  cylinders,  solenoid- 
operated  actuators,  and  a  control  panel  with  remote  enable  and  function  buttons.  The 
cylinders  had  a  volumetric  capacity  of  0.15  m®  arid  vented  through  a  top-mounted  2 b-cni 
(10-in)  butterfly  valve.  Trial  preparation  began  with  the  filling  of  one  or  both  cylinders  with 
SFe  gas  to  a  pressure  of  approximately  1 50  psig.  A  cylinder  pressure  of  1 50  psig  was  high 
enough  to  quickly  expel  the  contained  gas,  but  low  enough  to  prevent  the  SFe  frorn 
liquifying,  thereby  producing  a  puff  while  minimizing  unwanted  momentum  and  dense  gas 
effects  on  initial  puff  dimensipns.  Remote  control  buttons  vyere  used  to  enable  and  operate 
the  valves  either  together  or  separately.  Each  cylinder  was  instrumented  with  a  surface- 
mounted  thermometer  and  internal  cylinder  gage  pressure  transducer  to  provide  the 
measurements  required  for  calculation  of  the  released  mass.  The  release  system  was 
mounted  on  the  bed  of  a  5-ton  flatbed  truck.  This  truck  was  driven  to  one  of  the  desig¬ 
nated  release  locations  prior  to  the  beginning  of  each  trial  fill  procedure.  The  rdease  mass 
calculation  procedure  is  described  in  Section  3.1.3,  with  further  details  provided  in 
Appendix  E.  Table  2  in  Section  3  lists  the  mass  of  SF©  calculated  for  each  release. 
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RFCTION  3.  TRIAL  DATA  SUMMARIES 


3.1  DISSEMINATION  SUMMARIES 

3.1.1  Trial  Name  And  Time  Convention 

Dipole  Pride  26  trials  are  named  DSWAXX,  where  XX  is  a  sequential  numbering  of 
trials  (01  through  17).  All  SFg  gas  concentration  data  are  identified  by  trial  name.  Several 
of  the  trials  included  two  sets  of  disseminations,  the  second  release  being  staggered  90 
min  behind  the  first.  Consequently,  trial  names  are  supplemented  with  a  release  time 
convention  of  JJJhhmm.  The  release  time  convention  for  DP26  consists  of  the  three-digit 
Julian  date  (JJJ)  followed  by  the  hour  and  minute  of  the  dissemination  in  Pacific  Standard 
Time  (PST).  Dissemination  dates,  times,  and  locations  are  presented  in  Table  2. 

Two  time  conventions  are  used  with  DP26  meteorological  data.  The  MEDA  station 
data  collected  during  the  test  program  and  archived  on  CD  are  presented  in  PST,  with  the 
time  stamp  taken  from  the  end  of  the  averaging  period.  The  radiosonde  and  pibal  launch 
times  are  also  presented  in  PST.  In  contrast,  the  15-min  averaged  sonic  anemometer  data 
documented  in  this  report  and  on  CD  take  their  time  stamp  from  the  beginning  of  the 
averaging  period  and  are  presented  in  Universal  Coordinated  Time  (UTC),  which  is  8  hours 
ahead  of  PST.  For  example,  MEDA  station  data  for  the  period  1500-1515  PST  are  labeled 
1515,  while  sonic  data  from  the  same  period  are  labeled  2300  UTC. 

3.1.2  Dissemination  Procedures 

Trial  preparation  began  with  a  decision  to  disseminate  from  one  of  the  north  (N2, 

N3)  or  south  (S2,  S3)  positions.  The  truck-mounted  disseminator  system  was  then  driven 
to  that  position.  The  SFg  dissemination  cylinder  fill  procedure  began  upon  receipt  of  range 
safety  clearance  from  the  DOE  Nevada  Operations  Office  representative.  The  cylinder(s)  to 
be  used  on  a  trial  were  filled  to  a  nominal  gage  pressure  of  1 50  psi.  Cylinder  pressures  and 
temperatures  were  reported  when  the  fill  procedure  was  completed.  Small  SFe  leaks 
caused  by  temperature-induced  flexing  of  the  packing  material  between  the  cylinders  and 
valve  flanges  often  occurred  during  the  fill  procedure.  Leaks  were  minimized  by  tightening 
the  bolts  securing  the  valve  and  actuator  components  to  the  cylinder,  but  elevated  SFg 
background  concentrations  due  to  disseminator  leaks  are  apparent  in  the  sampler  data  for 
some  trials.  Table  2  includes  comments  concerning  trials  where  leaks  were  most  notice¬ 
able. 

3.1.3  Dissemination  Mass  Calculations 

Given  a  known  cylinder  volume,  the  cylinder  temperature,  and  internal  pressure,  it  is 
possible  to  calculate  the  mass  of  material  released.  Dr.  William  Espander  performed  a 
series  of  mass  calculations  based  on  the  Law  of  Corresponding  States,  a  Virial  equation, 
and  a  Martin-Hou  equation  (Mears  et  al,  1969).  While  the  results  from  these  three 
methods  were  similar,  Dr.  Espander  recommends  the  Martin-Hou  method  because  it  is 
based  on  actual  SFg  experimental  data.  Table  2  includes  disseminated  mass  calculated 
using  the  Martin-Hou  method.  Details  of  Dr.  Espander’s  DP26  calculation  procedures  are 
presented  in  Appendix  C. 
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3.1.4  Release  Dimensions 


An  IR  imager  positioned  100  m  east  or  west  (crosswind)  of  the  disseminator  pro¬ 
vided  puff  alongwind  and  vertical  dimension  information.  Release  duration  (usually  1-2  s)  is 
the  time  between  the  opening  of  the  first  valve  to  the  evacuation  of  the  last  cylinder.  The 
release  durations  were  subjectively  determined  from  visual  inspection  of  the  puff  release 
imagery.  Useable  puff  images  were  not  obtained  for  each  release,  and  those  that  were 
available  are  subject  to  uncertainties  in  interpretation.  Consequently,  this  report  provides 
no  trial-by-trial  source  dimension  information.  However,  there  was  enough  useful  data  to 
define  general  source  dimension  characteristics.  Momentum  of  the  exiting  gas  typically 
carried  the  puff  centroid  to  a  height  of  6  ±  2  m,  creating  SFe  puff  with  a  vertical  dimension 
of  4  ±  0.5  m.  The  initial  alongwind  puff  diameter  averaged  7.5  ±  2m.  Because  the  two 
cylinders  were  separated  by  less  than  2  m,  the  number  of  cylinders  used  during  a  release 
had  a  negligible  effect  on  source  dimensions. 

3.2  TRIAL  CONCENTRATION  SUMMARIES 

3.2.1  Detector  Placement  And  Sampling 

A.  Whole  Air  Samplers.  Thirty  whole  air  samplers  were  mounted  on  posts  1 .5  m  AGL 
at  nominal  300-m  intervals  along  each  sampling  line.  The  operators  of  the  six  van-mounted 
continuous  analyzers  also  performed  the  whole  air  sampler  programming,  two  being 
responsible  for  each  sampling  line.  Upon  notification  of  a  trial  start  time  by  the  test 
director,  operators  downloaded  to  each  whole  air  sampler  a  start  time  and  a  1 5-min 
sampling  interval.  The  two  sampling  lines  closest  to  the  release  point  were  programmed  to 
begin  sampling  at  the  projected  release  time  (actual  release  times  were  occasionally 
delayed),  while  the  furthest  line  of  samplers  was  programmed  with  a  30-min  delay.  Once 
the  sampling  began,  the  program  stepped  through  a  procedure  that  sequentially  filled  Bags 

1  through  1 2.  Sampling  for  each  trial  was  complete  upon  the  sealing  of  Bag  1 2  on  the 
sampling  line  furthest  from  the  source.  The  sealed  bags  were  taken  to  the  gas  laboratory 
for  analysis  at  the  end  of  each  trial. 

B.  Continuous  Analyzers.  Upon  completion  of  sampling  line  programming,  the  TGA- 
4000  operators  drove  their  vans  to  their  sampling  stations.  On  the  first  trial  (DSWA01 ), 
each  pair  of  vans  was  stationed  on  its  respective  sampling  lines  (Vans  3  and  4  on  Line  1 , 
Vans  1  and  6  on  Line  2,  and  Vans  2  and  5  on  Line  3).  However,  for  the  remainder  of  the 
trials,  a  decision  was  made  to  maximize  coverage  on  the  sampling  line  along  Pahute  Mesa 
Road  by  stationing  the  vans  as  indicated  in  Table  3.  This  arrangement  provided  real-time 
gas  concentration  readings  at  1500-m  intervals  across  Line  2  (see  Figure  1). 


Table  3.  Continuous  Analyzer  Van  Locations  by  Sampler  Station  Number. 

Van  1  Van  2  Van  3  Van  4  Van  5  Van  6 

230  224  212  218  206  201 


3.2.2  Sampler  Summaries 

Concentration  measurements  obtained  by  the  samplers  are  archived  on  a  CD  pro¬ 
duced  by  ARLFRD.  The  CD  contains  three  directories:  (1)  visualization  images  (Directory: 
IMAGES)  that  include  Yucca  Flat  terrain  overlaid  with  MEDA  station  winds  and  bargraph 
depictions  of  the  time  variations  of  SFg  concentrations;  (2)  a  directory  of  the  time-resolved 
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TGA-4000  data;  (3)  a  directory  of  the  whole  air  sampler  data  in  time-position  arrays.  The 
IMAGES  directory  presents  an  animated  depiction  of  the  puff  concentration  field  moving  in 
stages  through  the  sampling  lines.  The  other  directories  present  SFe  concentration  data  by 
position  and  bag  number  (whole  air  sampler  data)  or  by  van  number  and  time  (TGA-4000 
data). 


Whole  air  sampler  results  were  analyzed  to  obtain  normalized  concentration  maxima 
and  cloud  widths  (Oy)  for  each  trial  where  most  of  the  puff  crossed  the  sampling  lines. 
Normalization  was  accomplished  by  dividing  the  1 5-min  averaged  concentration  maxima  by 
source  strength  as  described  below.  The  results  are  presented  in  Table  4.  The  1 5-min 
time  resolution  available  from  the  whole  air  samplers  is  insufficient  to  resolve  alongwind 
cloud  growth  (a,)  dimensions  or  centroid  passage  time.  Continuous  analyzer  data  and 
remote  imagery  were  included  in  the  Dipole  Pride  26  test  to  provide  better  temporal 
resolution  of  puff  centroid  locations.  However,  because  the  continuous  analyzers  were 
stationed  only  along  Line  2  and  the  puff  imagery  has  not  been  analyzed,  alongwind  cloud 
dimensions  are  only  available  from  Line  2.  The  available  a,  results  are  summarized  in  Table 
5. 

3.2.3  Sampler  Quality  Control 

Watson  et  al.  (1998)  use  limit  of  detection  (LOD)  and  limit  of  quantitation  (LOQ)  to 
evaluate  TGA-4000  system  performance.  LOD  is  the  lowest  concentration  at  which  SFe 
can  be  detected,  defined  as  three  times  the  standard  deviation  of  measurements  at  zero 
concentration,  as  determined  by  analysis  of  signal  noise.  The  LOQ  is  defined  as  the 
minimum  concentration  measured  with  a  relative  error  within  ±  30  percent  at  the  95 
percent  confidence  level  (Taylor,  1987).  Watson  et  al.  (1998)  report  LODs  ranging  from  15 
to  42  pptv  and  LOQs  ranging  from  50  to  140  pptv  for  the  six  TGA-4000  units  deployed 
during  Dipole  Pride  26. 

Whole  air  sampler  quality  control  included  a  well-defined  series  of  procedural  steps 
designed  to  minimize  loss,  contamination,  and  mislabeling  of  samples  as  well  as  the  use  of 
duplicate  samples,  blanks,  and  spikes  to  define  accuracy  and  precision.  Sampler  quality 
control  and  performance  are  described  in  detail  by  Watson  et  al.  (1998).  A  total  of  19,688 
samples  were  collected.  Of  these,  1,386  were  rendered  unusable,  yielding  a  data  recovery 
rate  of  93  percent.  In  addition,  the  released  puff  apparently  missed  the  sampling  lines 
during  Trials  1,  2,  8,  and  10.  Consequently,  sampler  data  from  these  trials  were  not 
processed  and  are  not  included  in  the  trials  tabular  data.  This  reduces  the  available 
sampler  data  to  77  percent. 

LOD  baseline,  accuracy,  precision,  and  threshold  are  figures  of  merit  for  whole  air 
sampler  measurements.  Whole  air  sampler  baseline  and  accuracy  were  determined  using 
blanks  and  spikes.  A  blank  is  a  sample  collected  using  the  standard  sampling  protocol,  but 
with  the  bags  within  the  sampler  cartridge  containing  only  ultra  high  purity  air.  A  spike  is  a 
set  of  sampler  cartridges  containing  a  known  concentration  of  SFg  that  is  subject  to  the 
same  processing  as  test  data.  Analysis  of  455  blank  samples  yielded  a  mean  of  4  pptv 
with  a  standard  deviation  of  7  pptv.  This  result  produces  an  LOD  baseline  of  21  pptv. 
Spike  analysis  produced  a  mean  spike-to-standard  difference  of  6  pptv  with  a  standard 
deviation  of  15  pptv.  Watson  et  al.  (1998)  define  system  accuracy  at  ±  15  pptv.  Whole 
air  sampler  precision  was  determined  by  comparing  the  results  from  sets  of  duplicate 
samplers  positioned  on  sampling  Stations  115,  215,  and  315.  These  duplicate  samples 
were  divided  into  three  ranges  for  analysis:  0  to  100,  100  to  500,  and  greater  than  500 
pptv.  Watson  et  al.  (1998)  estimate  precision  at  the  95%  confidence  level  as  ±42%  for 
the  0  to  100  range,  ±32%  for  the  100  to  500  range,  and  ±1 2%  for  concentrations  greater 
than  500  pptv.  Threshold  considerations  include  the  effects  of  analytical  uncertainties  and 
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any  atmospheric  background  accumulation  of  SFg.  Accordingly,  samples  of  measurements 
taken  from  bags  exposed  to  ambient  air  not  in  the  vicinity  of  the  puffs  were  selected  for 
statistical  analysis.  These  samples  produced  a  mean  concentration  of  4.5  pptv  with  a 
standard  deviation  of  1.73.  There  was  also  no  noticeable  increase  in  the  SFg  background 
levels  from  the  beginning  to  the  end  of  the  trials  program.  The  background  threshold,  the 
lowest  likely  non-zero  reading  to  be  obtained  from  whole  air  samplers  operating  in  back¬ 
ground  air,  is  taken  as  three  standard  deviations  beyond  the  sample  rnean,  or  9.7  pptv. 

The  range  between  threshold  and  LOD  baseline  concentration  levels  (“10  to  21  pptv)  is  a 
"grey  area"  in  which  measurement  results  are  greater  than  those  expected  from  the  SFg 
background,  but  less  than  the  baseline  concentration  expected  at  the  edge  of  a  well- 
defined  SFe  puff.  Concentration  measurements  in  this  range  are  associated  either  with 
fugitive  emissions  from  a  leaking  disseminator,  or  with  the  periphery  of  a  very  diffuse  puff. 

3.2.4  Puff  Width  Estimates 

The  whole  air  sampler  data  available  for  each  trial  consist  of  samples  from  30 
crosswind  sampler  positions  along  each  of  the  three  lines.  Whole  air  samplers  mounted  at 
the  sampling  positions  provided  1 2  (Bags  1  through  1 2)  time-sequenced  concentration 
measurements,  generating  a  30  x  1 2  array  of  concentration  measurements  at  each 
sampling  line.  Figure  2  depicts  this  array  for  the  middle  sampling  line  of  Trial  DSWA1 2. 

The  1 5-min  sampling  intervals  set  for  each  bag  provided  insufficient  temporal  resolution  to 
estimate  centroid  passage  time  or  o,.  However,  the  300-m  sampler  line  spacing  usually 
provided  a  sufficient  number  of  above-threshold  concentration  measurements  along  each 
sampling  line  (a  minimum  of  six  is  desirable)  to  generate  histograms  of  the  lateral  puff 
concentration  distribution.  Gaussian  fits  to  these  histograms  were  used  to  determine  puff 
width  sigmas  (Oy).  "Best  estimates"  of  Oy  were  obtained  using  the  following  procedure: 

1 .  Define  a  coordinate  system.  The  DP26  field  test  domain  was  bounded  in  latitude 
by  the  north  and  south  puff  release  positions  (N2,  N3,  S2,  and  S3),  and  in  longitude  by 
sampling  positions  that  define  the  east  and  west  ends  of  each  sampling  lirie.  Watson  et  al. 
(1998)  provide  latitude  and  longitude  measurements  for  each  sampler  station  and  puff 
release  position.  The  position  of  a  puff  crossing  a  sampling  line  was  determined  by  the 
location  of  its  centroid  relative  to  the  sampler  positions,  where  the  centroid  was  defined  as 
the  center  of  mass  of  the  concentration  distribution.  To  determine  the  centroid  location,  a 
coordinate  system  was  established  along  each  line  with  a  local  origin  at  each  line  s  eastern¬ 
most  sampling  position  (130  for  Line  1,  230  for  Line  2,  330  for  Line  3).  All  other  samplers 
along  each  line  were  assigned  locations  at  nominal  300-m  increments  west  of  these 
origins.  This  arrangement  defined  a  lateral  distance  in  meters  from  the  local  origin  for 
every  sampler  position  along  each  line,  and  these  distances  were  used  with  the  sampler 
concentration  measurements  to  create  the  distance-concentration  histograms  shown  in 
Appendix  B.  The  histogram  centroid  distances,  determined  using  procedures  described 
below,  were  then  translated  back  to  latitude/  longitude  positions.  This  centroid  position 
information  was  used  with  release  position  information  to  determine  straight-line  distances 
from  the  puff's  source  location  to  its  position  as  it  crossed  each  sampling  line.  These 
distances  are  listed  under  "Distance  from  Source"  in  Table  4. 

2,  Determine  the  puff  histogram  concentration  maximum  (CJ.  For  each  30  X  12 
whole  air  sampler  concentration  array,  the  near-surface  puff  concentration  maxirnum  C„ 
was  identified  as  the  sampler  position  and  bag  number  near  the  array  concentration 
centroid  reporting  the  maximum  concentration.  Occasionally,  two  adjacent  bags  provided 
high  concentrations  (within  10  percent  of  each  other).  In  this  case,  the  true  likely 
passed  between  or  over  both  array  positions.  In  these  cases,  the  Oy  computation  procedure 
was  performed  using  each  set  of  bag  concentration  data.  Table  4  lists  the  1 5-min 
averaged  puff  C„  presented  in  pptv  and  as  normalized  (divided)  by  the  quantity  of  SFe 
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released  (see  Table  3)  using  a  conversion  procedure  described  by  ASTM  Standard  Practice 
D1 91 4-95  (ASTM,  1997B).  The  normalized  are  given  in  units  of  m'^  x  10’®. 

3.  Define  a  puff  vyidth  histogram.  The  bag  number  associated  with  the  puff 
concentration  maximurn  identifies  a  row  of  30  measurements  obtained  along  a  particular 
sampling  line  taken  during  the  same  1 5-min  period  as  the  measured  concentration  maxi¬ 
mum.  Concentration  measurements  and  position  information  from  the  selected  row  were 
used  to  construct  the  puff  width  histogram.  All  sampler  position  measurements  reporting 
above-baseline  concentrations  (in  excess  of  21  pptv)  within  this  time  period  were  used  in 
the  histogram.  Figure  2  illustrates  a  whole  air  sampler  concentration  array  obtained  during 
Trial  DSWA12,  Line  2  with  the  concentrations  used  to  define  puff  width  during  centroid 
passage  highlighted  in  bold  type.  A  total  of  40  useful  histograms,  identified  in  Table  4  by 
trial  name,  release  time,  sampling  line  number,  and  bag  number,  were  obtained  from  the 
Dipole  Pride  26  data  set.  Appendix  B  contains  additional  puff  width  histogram  information, 
including  conversions  from  bag  numbers  to  sampling  times. 

For  some  trials,  above-baseline  concentrations  were  measured  at  the  extreme  east 
and/or  west  sampler  positions,  indicating  that  the  puff  lateral  dimensions  likely  extended 
beyond  the  sampler  line.  If  the  missing  data  were  well  within  the  tail  of  the  distribution 
(greater  than  one  standard  deviation  from  the  centroid),  the  missing  part  of  the  histogram 
vvas  filled  in  using  estimates  of  centroid  shape  and  area  to  permit  statistical  analysis  of  the 
histogram.  Histogram  statistics  affected  by  this  procedure  are  enclosed  by  parentheses  in 
Table  4. 

4.  Plot  the  histogram  and  fit  a  Gaussian  curve.  The  crosswind  spread  of  diffusing 
material  is  traditionally  described  in  terms  of  a  Gaussian  Oy.  Although  material  released  into 
the  atrnosphere  initially  has  strong  concentration  gradients  at  the  boundaries,  entrainment  of 
clean  air  into  the  puff  and  diffusive  mixing  quickly  destroys  these  gradients  and  often  creates 
a  crosswind  concentration  field  that  resembles  a  Gaussian  distribution.  How  well  a  sampled 
concentration  field,  as  represented  by  its  concentration  histogram,  fits  a  Gaussian  distribution 
was  described  using  a  goodness-of-fit  parameter  and  the  histogram's  skewness  and  kurtosis. 

The  apparent  shape  of  a  puff's  crosswind  concentration  distribution,  as  represented 
by  its  concentration  histogram,  is  influenced  by  sampler  spacing  resolution,  and  the  nature  of 
the  flow  into  which  it  is  released.  When  a  small  puff  passes  over  only  a  few  crosswind 
samplers,  it  typically  produces  a  leptokurtic  (exaggerated  peakedness)  crosswind 
concentration  histogram  as  a  consequence  of  poor  spatial  resolution.  This  effect  is  evident  in 
some  of  the  concentration  histograms  obtained  from  puffs  crossing  the  nearest  sampling  line. 
As  the  puff  continues  to  expand  while  moving  downwind,  it  progressively  crosses  more 
sampling  stations  on  a  crosswind  line  and  the  sampler  spacing-induced  leptokurtic  effect 
diminishes.  With  progressively  greater  mixing,  the  histogram  eventually  assumes  a 
platykurtic  (exaggerated  flatness)  profile.  Wind  shear  and  divergence  can  also  skew  the 
concentration  distribution  and/or  create  secondary  peaks,  causing  the  concentration  field  to 
resemble  a  skewed  and/or  bimodal  distribution.  Therefore,  the  concentration-distance 
histograms  are  characterized  in  Table  4  using  a  "best-fit"  Oy  supplemented  by  the  coefficients 
of  skewness  and  kurtosis  and  a  measure  of  the  mean  square  departure  from  the  idealized 
Gaussian  distribution. 

The  histogram  analysis  procedure  began  with  keying  position  and  concentration  data 
into  the  Jandel  Scientific  SigmaPlot®  graphing  program.  A  transform  "GAUSq.XFM"  was 
written  to  operate  on  these  data.  GAUSq.XFM  calculates  the  histogram  area  (A),  puff 
centroid  (center  of  mass)  position  along  the  sampling  line,  "best-fit"  Oy,  and  the  mean  square 
error  of  the  Gaussian  curve  fit  to  the  histogram  data.  The  best-fit  Oy  is  the  one  that 
simultaneously  matches  histogram  area  (retaining  mass  continuity)  and  minimizes  the  least- 
squares  error  in  the  Gaussian  fit  to  the  histogram  data. 
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Figure  2.  Sulfur  hexafluoride  concentrations  in  pptv  arrayed  by  bag  number  and  whole  air 
sampler  location  along  Line  2  during  Trial  DSWA1 2  following  puff  releases  at 
0900  and  1030  PST  on  15  November  1996.  Puff  concentration  data  included  in 
puff  width  histograms  (described  in  3.2.4)  are  indicated  in  bold  type.  Bag  1 
provides  the  time-mean  concentration  during  the  first  1 5-min  time  period  after 
dissemination,  followed  successively  by  Bags  2  through  1 2. 


Table  4.  Dipole  Pride  26  Lateral  Dispersion  Summary  (time  in  PST). _ ; _ _ 

Trial  Release  Line  Bag  Distance  Observed  Maximum  Best  Fit  Coefficient  of  Normal  Crossing 

Name  date/time  Number  Number  from  Concentration  Oy  Skewness  Kurtosis  Departure  Angle 
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For  each  Gaussian  curve  of  area  A  and  dimension  Oy,  there  is  a  unique  concentration 
maximum  Cg  defined  by 


OyCg  =  A/(2p)°®  .  (3-1) 

Thus,  the  Oy  Cg  product  defines  a  family  of  curves  where  either  Sy  or  Cg  can  serve  as  an 
independent  variable.  If  the  histogram  C^  is  specified  as  Cg,  Equation  (3-1)  defines  a  Oyfor 
the  Gaussian  profile  that  matches  the  histogram  concentration  maximum.  However,  this  Sy 
may  not  provide  the  best  least-squares  fit  to  the  entire  histogram.  As  noted  above,  the  best 
least-squares  fit  was  obtained  by  iterative  convergence.  Because  the  SigmaPlot®  transform 
programming  does  not  provide  a  capability  to  converge  on  the  best  Gaussian  fit  to  histogram 
data,  the  best-fit  Cy  was  obtained  by  running  several  consecutive  solutions  for  GAUSq.XFM. 
GAUSq.XFM  includes  one  Oy  solution  with  the  histogram  input  as  Cg  (Solution  1 ),  a 
second  solution  which  calculates  a  Sy  using  the  standard  second-moment  statistical  metho¬ 
dology  (see  Kendall  and  Stuart,  1963),  and  a  third  "intermediate"  solution  based  on  a  mean 
Cg  derived  from  the  concentration  maxima  for  Solutions  1  and  2.  The  original  C^-based 
solution  usually  produced  a  relatively  small  Oy  for  the  Gaussian  curve  fitted  to  the  histogram 
maximum  concentration,  which  may  be  the  Cy  of  greatest  interest  if  predicting  the  maximum 
concentration  is  the  primary  modeling  objective.  The  second  (Kendall  and  Stewart)  solution 
typically  produced  a  larger  Oy  (and  smaller  Cg)  than  the  first  solution,  with  the  third  solution 
producing  Gy  and  Cg  intermediate  between  the  first  two.  Figure  3  illustrates  a  sample 
histogram  with  the  three  fitted  Gaussian  curves.  GAUSq.XFM  calculated  the  mean  square 
error  of  each  fit  to  the  histogram  and  normalized  it  by  the  total  histogram  concentration  SCj. 
Convergence  was  achieved  by  repetitively  solving  the  transform  using  Cg  estimates  that 
produced  progressively  smaller  mean  square  errors.  Sufficient  convergence  was  achieved 
when  the  difference  in  normalized  error  between  any  two  successive  GAUSq.XFM  solutions 
produced  error  differences  that  were  less  than  1 0  percent  of  the  minimum  normalized  error. 
For  each  puff.  Table  4  gives  the  best-fit  Gaussian  Oy  and  its  error  (listed  as  "normal  departure, 
%”).  An  additional  transform  (Gausfit.XFM)  provided  histogram  coefficients  of  skewness  and 
kurtosis  using  the  method  described  in  Kendall  and  Stuart  (1963).  These  statistics  are  also 
presented  in  Table  4.  Appendix  B  contains  printouts  of  the  Gausfit  and  GAUSq  transforms, 
plots  of  each  histogram,  and  the  best-fit  Gaussian  curve. 

Final  considerations  for  lateral  dimension  determination  include  detachment  of  the  puff 
centroid  from  the  surface  and  blending  of  the  puff  into  the  background.  When  the  bulk  of  the 
puff  is  detached  from  the  surface,  the  lateral  dimension  Oy  indicated  by  the  concentration 
measurements  along  the  sampling  lines  may  not  represent  the  actual  puff  Oy.  For  example, 
centroid  liftoff  was  observed  on  Line  3  during  release  3181400.  This  puff  lofted  the  first 
sampling  line,  producing  a  small  Cm  and  an  uncertain  Oy.  No  histograms  or  Oy  estimates  are 
presented  for  this  event  because  they  may  not  be  representative  of  values  at  the  elevation  of 
the  plume  centroid. 

t 

As  a  puff  proceeds  downrange,  it  becomes  progressively  more  dilute  and  eventually 
blends  into  the  background.  Accurate  determination  of  Oy  becomes  difficult  when  puff 
concentrations  diminish  to  magnitudes  between  threshold  and  baseline  concentration  levels. 
The  point  at  which  Oy  can  no  longer  be  determined  using  the  present  methodology  is  subjec¬ 
tive,  as  the  dilution  process  is  gradual  and  "grey  area"  puff  concentrations  affect  most  puff 
dimension  estimates  to  some  degree.  Because  puff  or  plume  borders  are  often  defined  at  10 
percent  of  the  maximum  concentration,  the  puff  Oy  estimates  in  this  report  are  based  on  only 
those  histograms  that  contained  a  maximum  concentration  in  excess  of  1 0  times  the 
threshold  (>97  pptv).  Other  puff  histogram  data  were  considered  too  dilute  for  accurate  Oy 
determination. 
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Trial  DSWA07,  Line  2.  Baa  3 


12  November  1996, 1330-1345  PST 

i  I  Concentration,  pptv 
— ^  Distance,  m  vs  s1  fit 
— ^  Distance,  m  vs  s2  fit 
— ^  Distance,  m  vs  s3  fit 


Figure  3.  Test  DSWA07,  Line  2,  Bag  3  (30-45  minutes  after  puff  release)  histogram  and 
Gaussian  curves  obtained  using  three  fit  methods  used  in  GAUSq.XFM. 
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5.  Compute  puff  crosswind  g,,  using  16-84  percentiles  from  the  histogram  CDF. 

The  sampler  position  and  concentration  information  used  to  calculate  the  best-fit  Sy  (Step  4 
above)  was  also  used  to  obtain  independent  histogram  Oy  estimates  based  on  the  1 and 
84*^  percentiles  of  the  concentration  cumulative  distribution  function  (CDF).  Because  68 
percent  of  the  area  under  a  normal  curve  falls  with  ±1  standard  deviation  of  the  mean,  the 
distance  between  the  locations  of  the  1 6*^  and  84*^  percentiles  spans  2  sigmas.  Dr.  Steven 
Hanna  (personal  communication)  considers  this  sigma  estimation  method,  which  requires 
interpolation  to  locate  the  1 S"*  and  84*^  percentile  positions  within  the  histogram,  to  be  more 
robust  than  the  moment  method.  Sigma  estimates  obtained  using  the  "Best-Fit"  and  "16-84" 
methods  were  nearly  identical  (within  ±10%)  for  most  cases,  and  differed  markedly  for  only 
two  cases  (DSWA06,  Line  3;  DSWA14,  Line  1)  where  the  concentration  distributions  were 
strongly  skewed  and  had  a  poor  fit  to  the  Gaussian  distribution. 

6.  Calculate  puff  travel  distances  and  crossing  angles  at  each  sampler  line.  The 
latitude  and  longitude  of  each  puff  centroid  crossing  position  was  calculated  using  the  puff 
centroid  line  crossing  position  determined  in  Step  2  above.  These,  along  with  the  latitude 
and  longitude  of  the  appropriate  release  position,  were  used  to  compute  the  distances  in 
meters  between  each  puff’s  origin  and  its  positions  as  it  crossed  the  sampling  lines.  The 
source-to-sampler  line  crossing  distances  are  shown  under  the  column  labeled  "distance 
from  source"  in  Table  4.  The  distances  were  originally  calculated  using  the  assumption  of 
a  spherical  earth.  Recalculation  with  corrections  for  the  oblate  spheroid  shape  of  the  earth 
produced  distances  diminished  by  a  factor  of  0.2  percent.  The  distances  presented  in 
Table  4  and  Table  5  include  this  correction. 

The  sampling  line  segments  at  the  point  of  puff  crossing  were  often  not  straight,  nor 
were  they  normal  to  the  wind  direction.  This  suggests  a  need  to  consider  the  apparent 
increase  in  a  puff’s  Oy  as  it  crosses  a  sampling  line  at  an  angle  of  less  than  90  degrees.  The 
puff  crossing  angle  was  determined  by  the  approach  angle  of  the  puff  centroid  to  the  line  that 
it  crossed.  Puff  crossing  angle  estimates  were  obtained  by  plotting  puff  centroid  tracks  from 
the  puff  origin  to  each  successive  sampling  line  crossing  position.  Crossing  angle  results 
presented  in  Table  4  only  estimate  the  puff  centroid  crossing  angle.  Most  of  the  crossing 
angles  were  found  to  be  large  (>45  degrees),  indicating  that  crossing  angle  effects  on  the 
magnitude  of  Oy  are  small.  The  Oy  information  provided  in  Table  4  are  not  adjusted  for 
crossing  angle. 

3.2.5  Puff  Alonowind  Dispersion  Estimates 

The  1 5-min  whole  air  sampler  time  resolution  was  insufficient  to  determine  puff  arrival 
and  passage  times  with  the  desired  accuracy.  However,  the  4-Hz  resolution  of  the  TGA- 
4000  measurements  provided  detailed  alongwind  histogram  data  that  permitted  a  precise 
determination  of  puff  centroid  times  of  arrival  and  transport  speeds.  The  difference  between 
the  time  of  a  puff’s  release  and  the  arrival  of  its  centroid  at  a  sampling  line  defines  puff  travel 
time,  and  the  straight-line  distance  traveled  by  the  puff  centroid  from  its  source  to  a  sampling 
line  divided  by  puff  travel  time  yielded  a  transport  speed.  Puff  release  time,  centroid  arrival 
time,  and  transport  speed  are  included  in  Table  5.  Puff  statistical  summary  data  are  missing 
for  Trial  DSWA1 1  because  the  gas  concentration  dilution  system  had  not  been  turned  on, 
causing  overranging  at  concentrations  in  excess  of  6934  pptv. 

Two  continuous  analyzers  were  initially  stationed  on  each  sampling  line,  but  after  Trial 
DSWA01  all  six  continuous  analyzers  were  stationed  at  positions  along  Line  2  to  maximize 
the  probability  of  the  puff  centroid  passing  near  at  least  one  of  these  instruments.  Thus, 
alongwind  puff  dimension  estimates  are  available  only  for  puffs  crossing  Line  2. 
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Table  5.  Dipole  Pride  26  Alongwind  Dispersion  Summary. 

Trial  Release  Position  Distance  PeaR  Best  Fit  Coetticient  ot  Normal  Centroid  Iransport 

Name  date/time  Number  from  Concentration  Sigma  t  Skewness  Kurtosis  Departure  Arrival  Speed 

(JJJhhmm)  Source  (ppt^)  (m'^XlO''^)  (s)  (ND)  (ND)  (%)  (hhmmss)  (m/s) 
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The  transforms  described  in  Section  3.2.4  and  a  methodology  similar  to  that  used  for  puff 
width  estimates  were  also  applied  to  alongwind  concentration  histograms.  The  coordinate 
system  used  to  define  puff  release  and  sampler  positions  was  also  used  to  define  the 
distance  frorn  the  source  to  centroid  arrival  at  Line  2.  The  best-fit  o,  defined  the  puff 
alongwind  dimension.  Table  5  presents  alongwind  puff  characteristics,  including  departure  of 
the  puff  time  histogram  data  from  the  best-fit  Gaussian  curve  and  histogram  skewness  and 
kurtosis. 

3.3  MICROMETEOROLOGICAL  SUMMARIES 
3.3.1  Statistical  Summaries 

Sonic  anemometer/thermometer  wind  and  temperature  component  measurements  at 
the  BJY  and  YFW  sites  were  reduced  to  1 5-min  averaged  statistical  summaries  for  the 
duration  of  each  trial.  This  averaging  produced  micrometeorological  statistics  for  periods  that 
correspond  to  the  MEDA  data  and  whole  air  sampler  averages.  The  statistical  summaries 
presented  In  Appendix  A  of  this  report  include  the  most  relevant  mean  and  second  moment 
quantities:  the  mean  wind  speed  and  direction,  horizontal  wind  angle  standard  deviation, 
alongwind  and  crosswind  velocity  variances,  and  fluxes  of  temperature  and  momentum. 
These  1 5-min  averaged  quantities  were  also  used  to  define  the  trial-averaged  stabilities  in 

Table  6. 


3.3.2  Roughness  Length  Estimates 

The  principal  application  of  surface  roughness  length  z,,  is  as  a  constant  in  the 
logarithmic  wind  profile  equation.  It  is  a  parameterization  of  the  degree  of  obstruction  to  free 
flow  presented  by  a  surface  to  air  moving  over  a  surface.  Both  the  size  and  spacing  of  flow 
obstructions  contribute  to  the  magnitude  of  z^.  When  flow  is  reasonably  steady  and  fully 
turbulent,  the  wind  speed  at  z  is  related  to  friction  velocity  u»  and  z/z^  by 

u,=  {ujk)mzlz„)  +  y/„,l  (3-1 

where  k  is  the  von  Karman  constant  (0.4  ±0.02)  and  y„,  is  a  stability-dependent  diabetic 
influence  function.  The  diabetic  influence  function  approaches  zero  as  heat  flux  approaches 
zero.  Thus,  the  logarithmic  wind  profile  equation  provides  a  means  of  solving  for  Zo  in  quasi- 
steady  near-neutral  conditions  when  rearranged  as  follows: 


=  z  e 


(-ku^  /  «,) 


(3-2) 


The  1 5-min  micrometeorological  summaries  obtained  from  the  sonics  at  BJY  and  YFW 
were  examined  for  cases  of  steady,  fully  developed  turbulent  flow.  Equation  (3-2)  was  used 
in  these  cases  to  solve  for  Zq  This  procedure  yielded  a  Zo  of  0.032  m  with  a  staridard 
deviation  of  0.012  m  at  BJY,  and  an  average  z^  of  0.045  m  with  a  standard  deviation  of 
0.026  m  at  YFW.  Because  YFW  is  on  the  edge  of  the  Yucca  Lake  salt  flat,  southerly  flow 
over  the  relatively  open  lake  bed  encounters  less  roughness  than  northerly  ^ow  over  the  more 
densely  vegetated  lake  bed  rim.  Northerly  flow  over  YFW  produced  a  Zq  of  0.06  Im  with  a 
standard  deviation  of  0.029  m,  while  southerly  flow  at  YFW  yielded  a  Zo  of  0.032  m  and  a 
standard  deviation  of  0.016  m. 


3.3.3  Boundary  Laver  Stability 


The  stability  of  the  atmospheric  layer  into  which  a  puff  is  released  can  be  used  to 
parameterize  the  turbulence  available  to  disperse  the  puff  material.  Stability  can  be 
parameterized  in  various  ways  ranging  from  a  simple  Pasquill-Gifford  methodology  (usually 
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based  on  wind  speed,  cloud  cover,  and  sun  angle)  to  the  Obukhov  length  L,  which  is 
inversely  proportional  to  the  ratio  of  buoyancy  and  shear  terms  in  the  turbulent  kinetic  energy 
equation.  Surface  weather  observations  and  MEDA  station  data  are  sufficient  for  the 
Pasquill-Gifford  methodology,  while  the  sonic  data  from  BJY  and  YFW  provide  the  fluxes  of 
heat  and  momentum  that  can  be  used  to  calculate  L.  Table  6  lists  the  wind  speed  and  cloud 
cover  information  required  for  the  Pasquill-Gifford  methodology  along  with  friction  velocities, 
temperature  fluxes,  and  estimates  of  L  derived  from  sonic  data. 

The  trial-averaged  friction  velocity,  temperature  flux,  and  Obukhov  length  estimates 
presented  in  Table  6  were  calculated  using  the  1 5-min  averaged  sonic  statistical  summaries 
provided  in  Appendix  A.  Although  the  1 5-min  average  is  somewhat  arbitrary  and  truncates 
some  of  the  low  frequency  motions  acting  on  the  puff,  this  averaging  period  is  expected  to 
capture  most  of  the  energy  involved  in  the  internal  mixing  of  puff  material.  The  scales  of 
motion  eliminated  by  this  "high  pass"  filtering  are  more  likely  involved  with  gross  puff 
movement  or  displacement. 

A  major  advantage  of  using  a  sonic  anemometer/thermometer  is  that  the  higher  order 
statistics  needed  to  characterize  the  state  of  the  atmospheric  boundary  layer  can  be  obtained 
directly  from  the  measurements  using  eddy  correlation  techniques.  However,  it  is  often 
difficult  to  obtain  statistically  stable  and  representative  momentum  and  temperature  fluxes, 
particularly  during  light  winds.  Busch  and  Panofsky  (1968)  recognize  this  problem  and 
consider  direct  momentum  flux  measurements  to  be  unreliable  when  the  friction  velocity  is 
less  than  0.32  m/s.  Fortunately,  momentum  flux  is  strongly  correlated  with  the  standard 
deviation  of  vertical  velocity  Biltoft  (1997)  finds  that  a  a„/u.  ratio  of  1 .3  can  provide 
reasonable  u.  estimates  under  low  wind  conditions.  This  ratio  was  used  to  estimate  u.  in  the 
low  wind  speed  cases  where  reliable  eddy  correlation  measurements  of  friction  velocity  were 
not  available. 

Temperature  flux  measurements  during  light  winds  are  also  problematic.  A  quiescent 
nocturnal  boundary  layer  is  characterized  by  small  vertical  motions  that  are  poorly  correlated 
with  small  temperature  fluctuations  punctuated  by  intermittent  turbulence  bursts.  A  majority 
of  the  vertical  heat  transfer  occurs  during  these  intermittent  bursts.  Although  the  long  term 
net  area  temperature  flux  must  be  downward,  individual  turbulence  bursts  observed  using 
measurements  at  a  single  site  may  produce  results  of  either  sign.  Consequent-ly,  there  can 
be  considerable  uncertainty  in  the  representativeness  of  tempera-ture  flux  measurements 
obtained  from  a  single  measurement  location  in  a  quiescent  nocturnal  boundary  layer.  These 
conditions  were  observed  at  YFW  during  the  early  morning  releases.  The  resultant 
temperature  flux  and  Obukhov  length  calculations  are  suspect,  especially  for  Trials  DSWA03 
and  DSWA06. 

3.4  MEDA  STATION  AND  UPPER  AIR  DATA 

MEDA  station  data  available  during  Dipole  Pride  26  include  1 5-min  measuremerits  of 
wind  speed,  wind  direction,  temperature,  and  pressure.  The  most  relevant  MEDA  stations 
(those  located  in  or  near  Yucca  Flat)  are  indicated  on  Figure  1  and  in  Table  1 .  Upper  air 
summaries  include  hourly  pibal  flights  interspersed  with  radiosonde  flights  taken  every  3 
hours.  The  pibal  flights  were  launched  from  DIGIPID  stations  at  BJY,  UCC,  and  CSE  and 
radiosonde  flights  were  launched  from  the  UCC  site.  MEDA  station,  pibal,  and  radiosonde 
flight  data  are  available  on  disk  and  on  CD. 

3.5  PUFF  IMAGERY 

Aerospace  Corporatiori  used  infrared  imagers  and  FTIRs  to  document  the  position  of 
each  puff  as  it  progressed  downwind  from  the  release  point.  In  some  cases,  the  puff  was 
tracked  well  beyond  the  most  distant  sampling  line.  Of  particular  relevance  to  this  data  set  is 
the  position  information  obtained  as  the  puff  centroids  crossed  the  sampling  lines.  This 
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information  is  needed  to  define  transport  wind  speeds  and  to  determine  whether  the  cenyoid 
was  in  contact  with  the  surface-based  sampler  arrays  or  passed  over  them.  Unfortunately, 
puff  imagery  was  not  available  before  the  completion  of  this  report. 
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SFCTION  4.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  Dipole  Pride  26  test  program  produced  a  comprehensive  set  of  long  range  puff 
dispersion  data  accompanied  by  extensive  meteorological  measurements.  Some  puffs  were 
tracked  to  distances  in  excess  of  20  km,  providing  an  opportunity  to  validate  both  the 
transport  and  diffusion  components  of  integrated  wind  field  and  dispersion  model  systems. 
This  report  presents  puff  Oy  and  a,  summaries  obtained  from  tracer  concentration 
measurements  along  with  trial  micrometeorological  statistics.  Data  sets  on  CD  ROM 
supplement  the  summaries  presented  in  this  report.  With  the  exception  of  unanalyzed  puff 
imagery,  this  report  and  the  associated  CD  ROM  constitute  a  fairly  complete  documentation 
of  the  Dipole  Pride  26  test  program.  Analysis  of  the  Aerospace  Corporation  puff  imagery 
could  provide  valuable  information  on  puff  vertical  dimensions  and  would  provide  independent 
estimates  of  puff  arrival  times  and  transport  speeds. 

The  1 5-min  averaged  micrometeorological  statistics  from  BJY  and  YFW  exhibit 
interesting  features  that  deserve  further  study.  The  fluxes  and  variances  sometimes  change 
by  more  than  an  order  of  magnitude  from  one  1 5-min  period  to  the  next  and  vary  greatly  ^ 
between  sites,  indicating  the  presence  of  terrain-induced  non-stationarity  and  inhomogeneity 
representative  of  dispersion  conditions  in  mountain-valley  desert  terrain.  Detailed  analysis  of 
micrometeorological  conditions  and  related  puff  behavior  are  beyond  the  scope  of  this  report, 
but  terrain-induced  differential  heating  and  cooling  appeared  to  exert  major  influences  on  puff 
behavior.  Differential  heating  of  mountain  slopes  created  local  flows  that  sometimes  drew 
puff  material  up  the  slopes  rather  than  along  the  major  axis  of  the  valley.  Also,  cold  air 
draining  from  the  north  often  pools  in  the  Yucca  Lake  basin  during  the  early  morning  hours. 
These  pools  of  confined  cold  air  stagnate  over  Yucca  Lake  until  post-sunrise  surface  heating 
destroys  them.  Puffs  released  into  this  early  morning  drainage  often  did  not  clear  Yucca  Lake 
until  several  hours  after  sunrise. 
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Table  6.  Dipole  Pride  26  Atmospheric  Stability  Summary.  _ 

Trial  Release  10  m  Wind  Speed  Cloud  Friction  Velocity  Temperature  Flux  Obukhov  Length 

Name  date/time  (m/s)  Cover  Height  (m/s)  (m  °K/s)  (m) 

(JJJhhmm)  BJY  YFW  (tenths)  (m)  BJY  YFW  BJY  YFW  BJY  YFW 

DSWA01  3091441  2J  O  9  7600  ;2467  .2766  !oTT3  -.0022  ^98  707 
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SECTION  5.  APPENDICES 


APPENDIX  A.  MICRQMETEOROLOGICAL  SUMMARIES 

The  micrometeorological  summaries  include  1 5-min  averages  of  sonic  anemometer/ 
thermometer  measurements  made  at  BJY  and  YEW.  The  measurement  height  was  1 0  m 
AGL.  The  summaries  include  average  wind  speed  and  direction,  horizontal  wind  angle 
standard  deviation  (oe),  alongwind  velocity  variance  ,  vertical  velocity  variance 
/  vertical  sonic  temperature  flux  (TTr),  and  momentum  flux  {u'w')  • 
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Trial:  DSWA17  Release  Julian  Date/Time:  3251330  Date/Time  (PST):  20  November  1330-1530  Source  Location:  S2 
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APPENDIX  B.  TRANSFORMS  AND  HISTOGRAM  PLOTS 


Appendix  B  contains  a  printout  of  the  GAUSq.XFM  and  Gausfit.XFM  transforms 
developed  using  the  Jandel  Scientific  Sigma  Plot®  program.  These  transforms  were  used 
to  produce  concentration  histograms  and  calculate  relevant  puff  dimension  statistics.  The 
transform  printouts  are  followed  by  plots  of  the  lateral  and  alongwind  puff  concentration 
histograms.  The  lateral  histogram  plots  include  a  best  fit  curve,  and  estimates  of  the  best 
fit  sigma  and  the  sigma  obtained  using  the  16-84  percentile  method  (see  Section  3.2.4). 
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jsv5R 

;GAUSq.XFM 

/Calculates  sigmas  for  histogram  data  and  provides  a 

/"best  fit”  sigma  from  a  range  of  three  sigma  estimates  calculated 

/for  Gaussian  fits  to  the  histogram  with  histogram  area  conserved. 

/  Calculations  are  done  on  user- 

/entered  class  mark  (position)  data  entered  in  column  1  and 
/ concentration  data  in  column  2 . 

/o  =  user-defined  arbitrary  origin 

/x=class  mark  (center  position)  of  each  histogram  column. 

/d=distance  from  each  x  to  o. 

/N=sum  of  concentrations  over  histogram. 

/y=concentration  (f)  for  each  histogram  column. 

/int=class  interval:  user-defined  histogram  width. 
min=  -2100  /user-defined  min  class  mark 
max=  5700  /user-defined  max  class  mark 
int  =  300  /  user-defined  class  interval 

o  =  300  /  user-defined  origin 

yO  =  113  /  user-entered  max  concentration  from  histogram 

pi=3. 1415926 

x_col=l 

y___col=2 

x=col(x_col)  /  user-entered  class  mark  in  column  1 

y=col  (y__col)  /user-entered  concentration  in  column  2 

N=total(y)  /sum  of  histogram  concentrations 

cell (3, 1) =”sum  cone:" 

cell(3,2)=  N 

cell (3, 3) ="hist  area:” 

A  =  N*int  /  histogram  area  calculation. 
cell(3,4)  =  A 

col  (4)  =  (x-o)/int  /  calculates  distance  d  from  user-selected  origin, 

col  (5)  =  col (4) *col (4)  /  calculates  d  squared. 

col (6)  =  col (2) *col (4)  /  product  of  d  and  concentration. 

col (7)  =  col (2) *col (5)  /  product  of  d**2  and  concentration. 

fd  =  total (col (6) ) 

/  histogram  mean  position,  xbar  calculation: 

xbar  =  o  +  int* (fd/N) 

cell (3,5)=  "xbar: " 

cell  (3,6)=  xbar  /  histogram  mean 

/first  (si)  histogram  standard  deviation  calculation. 

fd2  =  total (col (7) ) 

cell(3,7)=fd2 

si  =  int*sqrt ( (fd2/N) - ( (fd/N) **2) ) 
cell (3, 8)  =  "si,  1st  est” 
cell(3,9)  =  si 

/second  (s2)  histogram  sd  calculation. 

cell(3,10)="s2,2nd  est" 
s2=A/ (sqrt (2*pi) *y0) 
cell(3, ll)=s2 

/find  yl,  the  peak  concentration  from  the  data  set 
/that  produced  si  and  calculate  s3,  another  equivalent 
/area  sigma  between  si  and  s2 

yl=A/ (sqrt (2*pi) *sl) 
ynew=(y0+yl) /2 

s3=A/ (sqrt (2*pi) *ynew)  /s3  =  intermediate  sigma  estimate 

cell(3,12)=”s3" 

cell(3,13)=s3 

/Calculate  tau's,  ordinates,  and  fits  for  the  3  sigmas 
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col  {8)  =  {x-“xbar) /si  ;tau  calc  for  si 

col (9) = (x-xbar) /s2  ;tau  calc  for  s2 

col (10) = (x-xbar) /s3  ;tau  calc  for  s3 

col (11) =exp (- (col (8 ) **2/2) )  /ordinate  calc  for  si 

col (12) =exp (- (col (9) **2/2) )  /ordinate  calc  for  s2 

col (13) =exp (- (col (10) **2/2) )  /ordinate  calc  for  s3 

col(14)=col(ll)*int*sqrt( (N*(N-1) )/(2*pij )/sl  /fit  for  si 

col(15)=col(12)*int*sqrt ( (N*(N-1) )/(2*pi) )/s2  /fit  for  s2 

col(16)=col(13)*int*sqrt  (  (N*(N-1)  )./(2*pi)  )/s3  /fit  for  s3 

/calculate  total  least  square  differences: 

col (17)=sqrt ( (col (14) “col (2) ) **2)  /difference  fit  1  and  histogram 
col (18) =sqrt ( (col (15) -col (2) ) **2)  /difference  fit  2  and  histogram 
col (19) =sqrt ( (col (16) -col (2) ) **2)  /difference  fit  3  and  historgram 
lsl=total (col  (17) ) 
ls2=total (col  (18) ) 
ls3=total(col(19) ) 

lslp= (Isl/N) *100  /fit  1  least  sq  error  as  %  of  total  cone. 

Is2p= (Is2/N) *100  /fit  2  least  sq  error  as  %  of  total  cone. 

Is3p= (Is3/N) *100  /fit  3  least  sq  error  as  %  of  total  cone. 

cell(3,14)=”lsl  %:” 

cell (3^ 15) =lslp  /display  fit  1  Is  error  as  %  of  A 

cell(3, 16)="ls2  %:” 

cell (3, 17)=ls2p 

cell(3,18)=”ls3%:” 

cell(3,19)=ls3p 

/Find  the  best  fit  sigma: 

sig=if (lslp>ls2p, s2, si)  /select  the  min  of  Islp  and  ls2p 
sx={lslp, ls2p} 

1  smin=min  ( sx ) 

sigl=if (ls3p<lsmin, s3, sig)  /select  the  min  of  ls3p  and  sig 

cell(3,20)="SIGMA:” 

cell (3/ 21) =sigl 

/Computation  of  coefficients  of  skewness  and  kurtosis: 

/col (20)=col(5) *col (6) 

/ f d3=total ( col ( 20 ) ) 

/ col (21) -col (5) *col (7) 

/fd4=total(col(21) ) 

/m3=fd3-3*fd*fd2+2*fd**3  / skewness 

/m4=fd4-4*fd*fd3+6* ( (fd**2) *fd2-3*fd**4  / kurtosis 

/sk=m3/(sigl)**3 
/cell (3,22)-"Coeff .  SK:" 

/cell (3,23)=sk 
/ K=m4 / ( s igl ) *  *  4 
/cell(3,24)=”Coeff .  K:” 

/cell (3,25) -K 


B-3 


Gausfit .xfm 


jsv5R 

; GAUSFIT. XFM 

/Calculates  moments  and  coefficients  of  skewness  and  kurtosis  for  histo 
gram  data . 

/■Calculations  are  done  on  user- 

; entered  class  mark  (position)  data  entered  in  column  1  and 
/•concentration  data  in  column  2. 

;o  =  user-defined  arbitrary  origin 

;x=class  mark  (center  position)  of  each  histogram  column. 

;d=distance  from  each  x  to  o. 

;N=sum  of  concentrations  over  histogram. 

;y=concentration  (f)  for  each  histogram  column. 

;int=class  interval:  user-defined  histogram  width. 
min=13.9291  ;user-defined  min  class  mark 
rtiax=  14.0187  ;user-defined  max  class  mark 
int  =  .00006945  ;  user-defined  class  interval 

o  =  13.9  ;  user-defined  origin 

yO  =3083.9  ;  user-entered  max  concentration  from  histogram 

pi=3. 1415926 

x_col=l 

y_col=2 

x=col(x_col)  ;  user-entered  class  mark  in  column  1 
y=col(y_col)  /user-entered  concentration  in  column  2 

N=total(y)  ;sum  of  histogram  concentrations 

cell (3/ 1) ="sum  cone:" 

cell  (3,2)=  N 

cell (3, 3) ="hist  area:" 

A  =  N*int  ;  histogram  area  calculation. 
cell(3,4)  =  A 

col  (4)  —  (x— o) /int  ;  calculates  distance  d  from  user— selected  origin, 

col (5)  =  col (4 ) *col (4)  ;  calculates  d  squared. 

col (6)  =  col (2) *col (4)  ;  product  of  d  and  concentration. 

col (7)  =  col (2) *col (5)  ;  product  of  d**2  and  concentration. 

fd  =  total (col (6) ) 

ml=fd/N 

;  histogram  mean  position,  xbar  calculation: 

xbar  =  o  +  int*(fd/N) 

cell (3,7)=  "xbar:" 

cell  (3,8)=  xbar  ;  histogram  mean 

/•first  (si)  histogram  standard  deviation  calculation. 

fd2  =  total (col (7) ) 

m2=fd2/N 

/•  Computation  of  coefficients  of  skewness  and  kurtosis: 

col(8)=col(5)*col(6) 

fd3=total (col (8) ) 

m3=fd3/N 

col(9)=col(5)*col(7) 
f d4=total ( col ( 9 ) ) 
m4=fd4/N 
inm2=m2-ml 
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inm3=m3-3*ml*ra2+2*  (ml)  **3  /third  moment 

mm4=m4-4*ml*m3+6* ( (ml) **2) *m2-3* (ml) **4  ;4th  moment 

sk=mm3/ (mm2) **1 . 5  ;coeff  of  skewness 

cell(3, 9)="Coeff  SK:" 

cell(3,10)=sk 

K=mm4/ (mm2) **2 

cell(3,ll)="Coeff  K:" 

cell (3, 12)=K 

cell(3,13)=fd 

cell(3,14)=fd2 

cell(3,15)=fd3 

cell(3,16)=fd4 

cell (3, 17) =mm2 

cell (3, 18)=mm3 

cell (3, 19) =mm4 
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APPENDIX  C.  SULFUR  HEXAFLUORIDE  MASS  CALCULATION  PROCEDURE 

The  sulfur  hexafluoride  mass  calculation  procedure  described  in  Appendix  C  was 
contributed  by  Dr.  William  Espander  of  Logicon  RDA. 
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LOGICOM 


Logicon  RDA 

SYSTEMS  TECHNOLOGY  OPERATION 
2600  Yale  Boulevard  S.E. 

Post  Office  Box  9377 
Albuquerque,  NM  87119-9377 

Tel:  505  842-8911 
Fax:  505  242-2249 


6  January 

1997 

To: 

Dugway  Proving  Ground  (Dr.  C.  Biltoft,  WD-M) 
West  Desert  Test  Center  Meteorology  Division) 
Dugway,  UT  840222-5000 

From: 

William  R.  Espander 

Subject: 

Dipole  Pride  26  Mass  Calculation 

1.0  Summary 

A  quick  check  of  the  ideal  gas  assumption  used  to  calculate  the  mass  released  for  the  Dipole 
Pride  26  test  series’  was  made  using  the  law  of  corresponding  states.  This  check  indicated  a 
difference  in  the  calculated  mass  on  the  order  of  fifteen  percent  for  a  nominal  operating  point. 
A  more  in  depth  evaluation  was  conducted  to  verify  the  initial  calculation.  Three  independent 
methods,  the  law  of  corresponding  states,  a  Virial,  and  a  Martin-Hou^  equation  of  state  were 
used  that  yielded  similar  results.  Therefore,  it  is  suggested  that  a  form  of  an  equation  of  state 
other  than  ideal  gas  be  used.  I  recommend  that  the  Martin-Hou  equation  of  state  be  used  as 
it  is  anchored  to  experimental  data  in  the  range  of  interest  for  the  experiment.  The  results  of 
the  comparison  are  summarized  in  the  Table,  Figure  1,  and  Figure  2. 


2.0  Discussion 

The  validity  of  ideal  gas  assumption  used  to  calculate  the  mass  released  from  the 
dissemination  cylinder(s)  during  the  Dipole  Pride  26  test  series  was  evaluated  by  calculation 
of  the  compressibility  factor  for  a  ’nominal’  test  point. 


'  Chris  Biltoft,  FAX  to  Gary  Ganong  /  Bill  Espander,  19  Dec.  96,  "Dipole  Pride  26 
Preliminary  Mass  Calculations." 

^  W.  H.  Mears,  E.  Rosenthal,  and  J.  V.  Sinka,  "Physical  Properties  and  Virial  Coefficients 
of  Sulfur  Hexafluoride,"  J.  Phys.  Chem.,  vol  73„  pp  2254-2261,  July,  1969. 
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The  Law  of  Corresponding  States  holds  well  in  many  instances  for  molecules  that  are  not 
polar  or  hydrogen  bonded.  The  dipole  moment  for  sulfur  hexafluoride  is  zero^  indicating  the 
molecule  is  not  polar.  The  critical  temperature  and  pressure  for  sulfur  hexafluoride  are 
318.69  K  and  3.77  MPa,  respectively.  The  nominal  temperature  and  pressure  were  assumed 
to  be  300  K  and  1  MPa.  The  compressibility  factor  is  defined  as 

Z  =  +  (0  Z^*^(7’^,P^)  ,  where  is  the  spherical  molecule  term,  is  a 

deviation  function,  and  ©  is  the  Pitzer  acentric  factor.  The  value  for  the  compressibility 

factor  based  on  the  tables  and  constants  in  reference  2  is  0.874,  i.e.,  compressibility  should  be 
considered. 

The  van  der  Waals  equation  of  state  was  used  to  calculate  the  expelled  mass.  The  form  of 

van  der  Waals  equation  used  was  (p  +  n^alV  )  (V  -  nb  )  =  nRT  The  coefficients  are 
a  =  7.857  bar  1^  /  mole^  and  b  =  0.08786  d  /  mole. 


A  Virial  equation  of  state  may  be  defined  in  the  form  of  ^ _ L  .  The  first  virial 

R  T  ^  v‘ 

coefficient,  Bq,  is  unity  and  the  second  virial  coefficient,  Bj,  is  defined  as 

n 

P,  =  53  a.  [TJT  -  1  ]  ,  reference  4.  The  value  for  the  coefficients  are  given  as: 

i=l 

aj 

278.8 

646.8 
335.1 
-71.75 


1 
1 

2 

3 

4 


^  Reid,  Sherwood,  and  Prausnitz,  The  Properties  of  Gases  and  Liquids.  3^^  ed.,  McGraw-Hill 
Book  Co,  1977. 

*  David  R.  Lide,  CRC  Handbook  of  Chemistry  and  Physics.  72”*'  edition,  CRC  Press, 
1991  -  1992.  C-3 
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The  Martin-Hou  equation  of  state  is  a  variation  of  a  Virial  equation  of  state. 


p  ^  ^  A,  .  B,r  -  c,exp(-i:r/T.) 

M  (V  -  i>y 


,  T  is  in  K,  V  is  in  cc/g,  P  is  in  bars,  K  is 


6.88302200,  T,.  is  318.80  K,  and  b  is  0.32736730.  The  constants  are  given  as: 


i 

Ai 

Bi 

Q 

1 

0.0 

0.56926365 

0.0 

2 

-4.99043505  10^ 

0.54854082 

-2.37588665  10 

3 

4.12453944  10" 

-0.334003447 

2.81955047  10 

4 

-1.61292746  lO"' 

0.0 

0.0 

5 

-0.48996987 

0.109417750 

-3.08268133  10- 

The  mass  for  each  trial  was  calculated  using  the  four  equations  of  state.  The  results  of  these 
computations  are  summarized  in  the  table  and  displayed  graphically  in  the  figures. 

Figure  1  shows  the  calculated  mass  for  each  case.  A  case  is  defined  as  a  single  condition  for 
a  dispersion  cylinder.  A  trial  consists  of  one  or  two  cases,  depending  on  the  number  of 
cylinders  used.  Note  that  the  Virial  equation  of  states  give  a  fifteen  to  twenty  percent 
increase  in  mass  over  the  ideal  gas  law  or  van  der  Waals  equation  of  state. 

The  compressibility  factor  relates  the  equations  of  state  back  to  the  ideal  gas  law,  figure  2. 
The  Virial  equations  of  state  indicates  a  fifteen  to  twenty  percent  reduction  in  compressibility 
for  the  conditions  of  interest.  Case  1  shows  a  much  smaller  effect  because  the  pressure  is  a 
factor  of  ten  below  the  critical  point  compared  to  the  remainder  of  the  cases  that  are  a  factor 
of  three  below  the  critical  point. 

3.0  Conclusions 

A  cursory  evaluation  of  the  validity  of  using  the  ideal  gas  law  to  calculate  the  sulfur 
hexafluoride  mass  in  the  dispersion  cylinder  indicated  a  fifteen  to  twenty  percent  difference 
between  assumed  equation  of  states.  It  is  felt  that  this  difference  is  due  to  the  operating  point 
being  close  to  the  critical  point.  This  difference  was  verified  using  three  different  approaches, 
the  Law  of  Corresponding  States,  a  Virial  equation  of  state,  and  a  Martin-Hou  equation  of 
state.  The  values  calculated  using  a  van  der  Waals  equation  of  state  do  not  agree  with  the 
three  above  approaches.  Since  the  Martin-Hou  equation  of  state  is  based  on  experimental 
data  in  the  region  of  interest  for  the  Dipole  Pride  26  experiment,  it  is  suggested  that  this 
formulation  be  used  to  calculate  the  mass,  moles,  of  gas  released. 
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Trial-Name  Temperature  Pressure  Calculated-Mass  Total-Mass 

[K]  [Psig]  Ideal  van  der  Virial  Martin-Hou  Ideal  van  der  Virial  Martin-Hou 
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Dipole  Pride  26  Mass  Expelled  (concluded) 

Trial-Name  Temperature  Pressure  Calculated-Mass  Total-Mass 

[K]  [Psig]  Ideal  van  der  Virial  Martin-Hou  Ideal  van  der  Virial  Martin-Hou 


Dipole  Pride  26  Mass  Calculation 


<  5  > 


6  January  1997 


Os 

CO 

00 

CO 

Os 

os 

VO 

00 

O 

cs 

CN 

VO 

o 

p 

r-* 

CN 

CO 

CO 

CO 

1— < 

1—4 

d 

P 

d 

d 

d 

d 

<S 

cs 

cs 

cs 

cs 

VO 

CO 

1-H 

OS 

VO 

VO 

00 

r- 

00 

OS 

CO 

VO 

cs 

00 

CO 

wo 

(N 

t-H 

d 

d 

d 

d 

d 

d 

cs 

CN 

ca 

cs 

CN 

cs 

CO 

wo 

1-H 

CO 

CO 

VO 

WO 

VO 

VO 

CO 

CN| 

CO 

00 

os 

OS 

OS 

00 

00 

OS 

1—4 

1-H 

1-H 

1-H 

v-H 

1-H 

xf 

WO 

CO 

CO 

wo 

VO 

VO 

CO 

CNl 

p 

OS 

00 

Os 

Os 

OS 

00 

00 

VO 

ON 

OS 


00 

wo 

wo 

VO 

VO 

CO 

Os 

Os 

CO 

VO 

^H 

wo 

o 

CO 

cs 

^H 

Q 

wo 

wo 

VO 

wo 

o 

r- 

wo 

CM 

o 

00 

o 

00 

oo 

o 

o 

O) 

p 

oq 

OS 

OS 

OS 

OS 

00 

OS 

00 

00 

OS 

OS 

00 

00 

00 

os 

wo 

40 

VO 

VO 

o 

O 

VO 

^-4 

VO 

o 

CO 

(N 

1—4 

o 

WO 

wo 

00 

wo 

O 

wo 

T— 4 

CN 

o 

oq 

o 

oq 

oq 

p 

O 

p 

p 

oq 

OS 

Os 

OS 

Os 

00 

OS 

00 

00 

OS 

OS 

00 

00 

00 

cn 

CO 

CO 

o 

r- 

o 

1—4 

1-H 

o 

o 

o 

p 

d 

d 

vd 

d 

os 

d 

d 

OS 

d 

d 

d 

wo 

wo 

wo 

wo 

wo 

wo 

wo 

wo 

WO 

wo 

WO 

1-H 

1—4 

1^ 

1-H 

1-H 

1—4 

ir-4 

1-H 

1—4 

1-H 

1-H 

1—4 

1—4 

oq 

P 

p 

p 

1-H 

p 

1-H 

p 

p 

o 

CO 

p 

00 

vd 

1—4 

r4 

(> 

wo 

1—4 

ci 

vd 

s 

os 

Os 

os 

o 

o 

o 

o 

o 

o 

o 

O 

o 

O 

(M 

cs 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

o 

o 

o 

o 

o 

Q 

o 

o 

CO 

o 

o 

CO 

o 

CO 

CO 

1-H 

CO 

CM 

CO 

CM 

CO 

1—4 

^H 

1-H 

1-H 

1—4 

CO 

CO 

wo 

wo 

cs 

CM 

CM 

CM 

CM 

CM 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

C-6 


g:\clltrl\pride\dp26\dp26niass.ltr  6  January  1997 


Dipole  Pride  26  Mass  Calculation 


<  6  > 


6  January  1997 


Dipole  Pride  26  Mass  Calculation  <  7  >  6  January  1997 


Sincerely, 


cc:  LTC  A.  J.  Kuehn  (DSWAAVELE) 

G.  Ganong  (LRDA/ABQ) 

T.  Mazzola  (LRDA/TGV) 
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